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Despite the advances in the development of numerical methods analytical approaches play a key 
role on the way towards a deeper understanding of strongly interacting systems. In this re- 
gards, renormalization schemes for Hamiltonians represent an important new direction in the 
field. Among these renormalization schemes the projector-based renormalization method (PRM) 
reviewed here might be the approach with the widest range of possible applications: As demon- 
strated in this review, continuous unitary transformations, perturbation theory, non-perturbative 
phenomena, and quantum-phase transitions can be understood within the same theoretical frame- 
work. This review starts from the definition of an effective Hamiltonian by means of projection 
operators that allows the evaluation within perturbation theory as well as the formulation of a 
renormalization scheme. The developed approach is then applied to three different many-particle 
systems; At first, we study the electron-phonon problem to discuss several modifications of the 
method and to demonstrate how phase transitions can be described within the PRM. Secondly, to 
show that non-perturbative phenomena are accessible by the PRM, the periodic Anderson is in- 
vestigated to describe heavy-fermion behavior. Finally, we discuss the quantum-phase transition 
in the one-dimensional Holstein model of spinless fermions where both metallic and insulating 
phase are described within the same theoretical framework. 
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I. INTRODUCTION 

During the last three decades the investigation of phe- 
nomena related with strongly interacting electrons has 
developed to a central field of condensed matter physics. 
In this context, high-temperature superconductivity and 
heavy-fermion behavior are maybe the most important 
examples. It has been clearly turned out that such sys- 
tems require true many-body approaches that properly 
take into account the dominant strong electronic corre- 
lations. 

In the past, many powerful numerical methods like ex- 
act diagonalization numerical renormalization group 
(0), Quantum Monte-Carlo ([!), the density-matrix renor- 
malization group or the dynamical mean-field the- 
ory (^) have been developed to study strongly correlated 
electronic systems. In contrast, only very few analytical 
approaches are available to tackle such systems. In this 
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regard, renormalization schemes for Hamiltonians devel- 
oped in the nineties of the last century @; 0; H) represent 
an important new direction in the field where renormal- 
ization schemes are implemented in the Liouville space 
(that is built up by all operators of the Hilbert space). 
Thus, these approaches can be considered as further de- 
velopments of common renormalization group theory (9) 
that is based on a renormalization within the Hilbert 
space. 

In this review we want to discuss the projector-based 
renormalization method [PRM, Ref. that shares some 
basic concepts with the renormalization schemes for 
Hamiltonians mentioned above (0; 0; i). AU these ap- 
proaches including the PRM generate effective Hamilto- 
nians by applying a sequence of unitary transformations 
to the initial Hamiltonian of the physical system. How- 
ever, there is one distinct difference between these meth- 
ods: Both similarity renormalization (@; 0) and Wegner's 
flow equation method (8) start from a continuous for- 
mulation of the unitary transformation by means of a 
differential form. In contrast, the PRM is based on dis- 
crete transformations so that a direct link to perturbation 
theory can be provided. 

This review is organized as follows: 

In the next section we discuss the basic concepts of the 
PRM: We introduce projection operators in the Liouville 
space that allow the definition of an effective Hamilto- 
nian. If these ingredients are combined with unitary 
transformations one can derive a new kind of pertur- 
bation theory that is not restricted to the ground-state 
but also allows to investigate excitations. (To illustrate 
this point we briefly discuss the triplet dispersion rela- 
tion of a dimerized and frustrated spin chain in the Ap- 
pendix.) However, this perturbation theory is not the 
focus of this review and can be considered as an inter- 
esting side-product of the development of the PRM, a 
renormalization scheme based on the same ingredients. 
To illustrate the method in some detail, the exactly solv- 
able Fano- Anderson model is considered. 

Improving our previous publications on the PRM, we 
show here the relation of the PRM to Wegner's flow equa- 
tion method (8) for the first time. It turns out the latter 
method can be understood within the framework of the 
PRM by choosing a complementary unitary transforma- 
tions to generate the effective Hamiltonian. For demon- 
stration, the Fano-Anderson model is solved with this 
approach, too. 

As a more physical example, the electron-phonon in- 
teraction is studied in Sec. lIIII In particular, the PRM is 
compared in some detail with the flow equation method 
(0) and the similarity transformation |g; 0). Further- 
more, we introduce a possible modification of the PRM 
that allows to derive block-diagonal Hamiltonians, and 
we discuss in some detail the freedom in choosing the gen- 
erator of the unitary transformation the PRM is based 
on. Finally, we show how phase transitions can be stud- 
ied within the PRM by adding symmetry breaking fields 



to the Hamiltonian. 

In Sec. IIVI the PRM is applied to the periodic Ander- 
son model to describe heavy-fermion behavior. Whereas 
the famous slave-boson mean-field theory (fill : fl^ ob- 
tains an effectively free system consisting of two non- 
interacting fermionic quasi-particles, here the periodic 
Anderson model is mapped onto an effective model that 
still takes into account electronic correlations. Thus, in 
principle both mixed and integral valence solution can be 
found. However, here we restrict ourself to an analytical 
solution of the renormalization equations that is limited 
to the mixed valence case. 

As third application of the PRM the one-dimensional 
Holstein model of spinless fermions is discussed. It is well 
known that the system undergoes a quantum phase tran- 
sition from a metallic to a Peierls distorted state if the 
electron-phonon coupling exceeds a critical value. First, 
for the metallic state we discuss the crossover behavior 
between the adiabatic and anti-adiabatic case in Sec. [V] 
All physical properties are shown to strongly depend on 
the ratio of phonon and hopping energy in the system. 
In Sec. IVI[ a unified description of the quantum-phase 
transition is given for the one-dimensional model in the 
adiabatic case. 

Finally, as a second example for a quantum phase tran- 
sition, we discuss in Sec. IVIII the competition of charge 
ordering and superconductivity in the two-dimensional 
Holstein model. Based on the PRM both charge den- 
sity wave and superconductivity are studied within one 
theoretical framework. 

We summarize in Sec. I Villi 

II. PROJECTOR-BASED RENORMALIZATION 
METHOD (PRM) 

In this section we introduce the concepts of the PRM 
jio!) where we particularly pay attention to a general no- 
tation that is used throughout the review for all applica- 
tions of the approach. 

We define projection operators of the Liouville space 
and define an effective Hamiltonian where, in contrast to 
common approaches, excitations instead of states are in- 
tegrated out. In this way, not only a perturbation theory 
is derived but also and more important a renormaliza- 
tion scheme (that we call PRM in the following) is estab- 
lished which allows to diagonalize or at least to quasi- 
diagonalize many-particle Hamiltonians. As an illustra- 
tive example, the exactly solvable Fano-Anderson model 
is discussed. 

The PRM is based on a sequence of finite unitary 
transformations whereas Wegner's flow equations start 
from a continuous formulation of unitary transformations 
by means of a differential form. It turns out that such 
a continuous transformation can also be understood in 
the framework of the PRM if a complementary choice for 
the generator of the unitary transformation is used and 
infinitely small transformation steps are considered. To 
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discuss the differences between tfie two formulations of 
the PRM in more detail, we also solve the Fano- Anderson 
model using the developed continuous approach. 

A. Basic concepts 

The projector-based renormalization method (PRM) 
(fiol ) starts from the usual decomposition of a given many- 
particle Hamiltonian, 

"H = Ti-o + TLi , 

where the perturbation Tii should not contain any terms 
that commute with the unperturbed part Hq. Thus, the 
interaction Hi consists of the transitions between eigen- 
states of Ti.0 with corresponding non-zero transition en- 
ergies. The presence of Tii usually prevents an exact so- 
lution of the eigenvalue problem of the full Hamiltonian 
Ti. so that suited approximations are necessary. 

The aim is to construct an effective Hamiltonian Ti.x 
with a renormalized 'unperturbed' part Tin, a a-nd a re- 
maining 'perturbation' Ti.i \ 

Ti-x = Hq^x+Hi^x (2.1) 

with the following properties: 

(i) The eigenvalue problem of the renormalized Hamil- 
tonian Ho,x is diagonal 

■Ho,x\n^) = E^\n^) 
with A-dependent eigenvalues and eigenvectors 

(ii) The effective Hamiltonian Hx is constructed in such 
a way so that (measured with respect to Hq.x) all 
non-diagonal contributions with transition energies 
larger than some cutoff energy A vanish. 

(iii) Tix has the same eigenvalues as the original Hamil- 
tonian Ti.. 

The eigenvalue problem of Tig, a is crucial for the con- 
struction of Hx because it can be used to define projec- 
tion operators, 

Pa^ = Y,\n^){m^\{n^\A\m^) (2.2) 

xeix-\E::-E^\) 

Qa = 1-Pa. (2.3) 

Note that neither \n^) nor |m'^) need to be low- or 
high-energy eigenstates of Ho.x- Pa and Qa are super- 
operators acting on operators A of the Hilbert space of 
the system. Thus, Pa and Qa can be interpreted as pro- 
jection operators of the Liouville space that is built up by 
all operators of the Hilbert space. Pa projects on those 



parts of an operator A which only consist of transition 
operators \n^){m^\ with energy differences \E^ — E^ \ less 
than a given cutoff A, whereas Qa projects onto the high- 
energy transitions of A. 

In terms of the projection operators Pa and Qa the 
property of Hx to allow no transitions between the eigen- 
states of 7^0. A with energies larger than A reads 

Hx = PxHx or QaT^a = 0. (2.4) 

For an actual construction of the effective Hamiltonian 
we now assume that the effective Hamiltonian Hx can 
be obtained from the original Hamiltonian by a unitary 
transformation, 

Hx = e^^ H e-^^ . (2.5) 

which shall automatically guarantee that condition (iii) 
above is fulfilled. 

In the following the evaluation of the effective Hamil- 
tonian (j2.5p is done in two ways: At first a perturba- 
tive treatment is derived. After that we develop a much 
more sophisticated renormalization where we interprete 
the unitary transformation of Eq. (j2.5p as a sequence of 
small transformations. The projector-based perturbation 
theory discussed in the next subsection is important for 
the understanding of the renormalization scheme derived 
later. However, the main focus of this review is the PRM. 

B. Perturbation theory 

In the following we evaluate the effective Hamiltonian 
Hx in perturbation theory. For this purpose the effec- 
tive Hamiltonian Hx from Eqs. (|2.4p and (|2.5p is sim- 
plified in a crucial point: The projection operators are 
now defined with respect to the eigenvalue problem of 
the unperturbed part of the original Hamiltonian Hq, 

Ho\n) = i?Jn). 

Thus, these projection operators differ from the formerly 
defined projectors Pa and Qa and can be written as fol- 
lows 

Pa^ = ^\n){m\{n\A\m) (2.6) 

xe(A-|£;„ -£;„!), 

Qa = 1-Pa. (2.7) 

The renormalized Hamiltonian Hx is now obtained from 
the unitary transformation (|2.5p . 

Hx = Pa^a = H e-^^ , 

where Xx is the generator of this transformation. To find 
Xx, we employ the modified condition (|2.4p : All matrix 
elements of Hx for transitions with energies larger than 
A vanish, i.e. 

Qa^a = (2.8) 
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First we expand Tix with respect to X\ 



Ti-x — Ti. + [X\ , Ti] 



1 
2! 



[Xx,[Xx,H]] (2.9) 



+ -[Xx,[Xx,[Xx,n]]] + ... . 

and assume that the generator X\ can be written as a 
power series in the interaction Tii, 



X, 



X 



(1) 



X 



(2) 



X 



f + (2.10) 



Thus inserting f^JU\i in Eq. (1^ . the effective Hamilto- 
nian Ti^ can be rewritten as a power series in the inter- 
action Til 



Ux 



Ti-o + "Hi 



X^^KHo 



1 

2! 



(1) 

A ' 



(2.11) 



o{nl) 



The contributions to the generator of the uni- 

tary transformation can successively be determined by 
employing Eq. (|2.8p . One finds 



^xX'x' 



1 

Lo 



(Qa^i) , 



Qa4'^ = -— Qa 



1 

2Lo 



Y^Qa 



(QA^l),i-(QAWl) 
(PAWl),:j^(QA^l^ 



(2.12) 
(2.13) 



Here, Lq is the Liouville operator of the unperturbed 
Hamiltonian TYq which is defined by \jqA = \HqtA\ for 
any operator variable A. 

As one can see from (|2.12p and (|2.13[) . no information 
about the low-energy part P^Xx of the generator Xx can 
be deduced from ()2.8|) . Therefore, we set for simplicity 



'(1) _ 



p,x, = VxX\'> = p,xr^ =•••= 



A-^A 



(2) _ 



(2.14) 



Inserting Eqs. ((2T2)) . ((2T3)) . and ((2T4)) into the power 
series (|2.1ip for T-L\, the desired perturbation theory is 
found, 



'Hx = Tio + Pa^i 



1. 



(QaWi):1^(QaWi; 



(P,7^i),— (Q^Hi) 



0(7^?), (2.15) 



which can easily be extended to higher order terms. Note 
that the correct size dependence of the Hamiltonian is au- 
tomatically guaranteed by the commutators in Eq. (j2.15p . 
The limit A —> is of particular interest because in this 
case the complete interaction Hi is integrated out. 

Usual perturbation theory derives effective Hamiltoni- 
ans that are only valid for a certain range of the system's 
Hilbert space. In contrast, Tix-, as derived above, has no 
limitations with respect to the Hilbert space so that it 



can also be used to study excited states. To illustrate 
this important aspect of our projector-based perturba- 
tion theory, we discuss the dimerized and frustrated spin 
chain in the appendix. 

At this point we would like to note that Eq. (|2.15p can 

(2) 

also be derived in a different way. It turns out that X^ ' 
is only needed to fulfill the requirement Tlx = Pa^^ 
we restrict ourselves to second order perturbation theory. 

(2) 

Thus, in this case X^ ' can be set to if the projector 
P;^ is applied to the right hand side of Eq. (|2.1ip . 



Hx - 7^o + Pa^i+Pa 
+ P ' ' 



X^^\Tio 



(2.16) 



-^A ' '^l 



Ip 

2!^ 



X 



(1) 

A ' 



X^^\Ho 



It is easy to proof that Eq. (|2.16p again leads to the result 
Eq. if and (^1^ is used. 

In Appendix [XI the developed perturbation theory 
(|2.15p is applied to the dimerized and frustrated spin 
chain where ground-state energy and triplet dispersion 
relation have been calculated. 

A perturbation theory based on Wegner's flow equa- 
tions (S), that also allows a description of the complete 
Hilbert space, has been derived in Refs. ITsI andH. How- 
ever, this approach requires an equidistant spectrum of 
the unperturbed Hamiltonian Tip- In contrast, the per- 
turbation theory presented here can be applied to sys- 
tems with arbitrary Hilbert space, and has similarities to 
a cumulant approach to effective Hamiltonians l\15t \. 



C. Stepwise renormalization 

In the previous subsection the effective Hamiltonian 
Hx as defined by Eqs. (|2.4p and (|2.5p has been evaluated 
within a new kind of perturbation theory. However, if the 
unitary transformation (j2.5p is interpreted as a sequence 
of unitary transformation a renormalization scheme can 
be developed based on the same definition of the effective 
Hamiltonian. Because again the projection operators Pa 
and Qa play a key role we call the derived method (fioi ) 
projector-based renormalization method (PRM). 

Let us start from a renormalized Hamiltonian Hx = 
T^o,x + Hi,x that has been obtained after all transitions 
with energy differences larger than A have already been 
integrated out. Of course, Ho,x and Hi.x will differ from 
the original Hq and Hi- Furthermore, we assume Hx has 
the properties (i)-(iii) proposed in subsection III. Al 

Now we want to eliminate all excitations within the 
energy range between A and a smaller new energy cutoff 
A — AA. Thereby we use a unitary transformation. 



H 



(A-AA) 



Hx e 



(2.17) 
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so that the effective Hamiltonian Hx-ax has the same 
eigenspectrum as the Hamiltonian H\. Note that the 
generator X\,ax needs to be chosen anti-Hermitian, 

-''^A.AA = ^-'^A AA' ensure that Hx-ax is Hermitian 
when Tix was Hermitian before. To find an appropriate 
generator Xx^ax of the unitary transformation, we em- 
ploy the condition that Tlx has (with respect to TYq.a) 
only vanishing matrix elements for transitions with ener- 
gies larger than A, i.e. Qx^x — 0. Similarly, also 

Q(A-AA)^(A-AA) = (2-18) 

must be fulfilled, where Q(a-aa) is now defined with re- 
spect to the excitations of 7io,(A-AA)- 

In principle, there are two strategies to evaluate 
Eqs. (PT7)) and (PTT5)l : The first uses perturbation the- 
ory as derived in subsection IH.BI In this case Hfx-AX) 
can be written as 

W(A-AA) = (2.19) 
= ^0,A + P(A-AA)^1,A + P(A-AA) [-'^A,AA, Wo,a] 
+ P(A-AA) [-'^A.AA, Wi,a] 

+ 2^(X-AX) [Xx,AX, [Xx,AX,'Ho,x]] + C>{HI x)- 

The generator Xx ax has to be chosen corresponding to 

Eq. arm . 

Q(x-AX)Xx.Ax = [Q(a-aa)^i,a] H (2.20) 

i^o.x 

For details of the derivation we refer to subsection 
III.Bl This approach has been successfully applied to the 
electron-phonon interaction to describe superconductiv- 
ity iS). 

Alternatively, one can also start from an appropriate 
ansatz for the generator in order to calculate H(a_aa) 
in a non-perturbative manner (f2l|). An ansatz for the 
generator with the same operator structure as Eq. (|2.20p 
is often a very good choice. This approach has been ap- 
plied to the periodic Anderson model to describe heavy- 
fermion behavior (U; 12^ ) . 

It turns out that the second strategy has the great 
advantage to successfully prevent diverging renormaliza- 
tion contributions. However, in both cases, Eqs. (|2.17p 
and (|2.18p describe a renormalization step that lowers 
the energy cutoff of the effective Hamiltonian from A to 
A— AA. Consequently, difference equations for the Hamil- 
tonian Tix can be derived, and the resulting equations 
for the A dependence of the parameters of the Hamilto- 
nian are called renormalization equations. By starting 
from the original model Ti. —: Ti.x=A the Hamiltonian is 
renormalized by reducing the cutoff A in steps AA. The 
limit A ^ provides the desired effective Hamiltonian 
7ix=o ='■ without any interaction. Note that the re- 
sults strongly depend on the parameters of the original 
Hamiltonian Ti.. 



D. Generator of the unitary transformation and further 
approximations 

It turns out that the generator Xx^ax of the unitary 
transformation is not yet completely determined by Eqs. 
(I2.17P and (|2.18p . Instead, the low-energetic excitations 
included in Xx.ax, namely the part P(a-aa)-^a,aa, can 
be chosen arbitrarily. The result of the renormalization 
scheme should not depend on the particular choice of 
P(A-AA)-'^A,AA as long as all renormalization steps are 
performed without approximations. However, approxi- 
mations will be necessary for practically all interacting 
systems of interest so the choice P(a-aa)-'^a,aa becomes 
relevant. If P(a-aa)^a,aa = is chosen the minimal 
transformation is performed to match the requirement 
(|2.18p . Such an approach of " minimal" transformations 
avoid errors caused by approximations necessary for ev- 
ery renormalization step as much as possible. Note that 
in order to derive the expression (|2.19p this choice of 
P(A-AA)-^A,AA was used. However, in particular cases a 
non-zero choice for P(a-aa)-^a,aa might help to circum- 
vent problems in the evaluation of the renormalization 
equations. 

In general, new interaction terms can be generated in 
every renormalization step. This might allow the investi- 
gation of competing interactions which naturally emerge 
within the renormalization procedure. However, actual 
calculations require a closed set of renormalization equa- 
tions. Thus, often a factorization approximation has to 
be performed in order to trace back complicated op- 
erators to terms already appearing in the renormaliza- 
tion ansatz. Consequently, derived effective Hamiltoni- 
ans might be limited in their possible applications if im- 
portant operators have not been appropriately included 
in the renormalization scheme. 

If a factorization approximation needs to be performed 
the obtained renormalization equations will contain ex- 
pectation values that must be calculated separately. In 
principle, these expectation values are defined with re- 
spect to Tix because the factorization approximation was 
employed for the renormalization step that transformed 
Hx to 7i(A-AA)- However, Hx still contains interactions 
that prevent a straight evaluation of required expecta- 
tion values. The easiest way to circumvent this difficulty 
is to neglect the interactions and to use the diagonal un- 
perturbed part 7^0, A instead of Hx for the calculation of 
the expectation values. This approach has been success- 
fully applied to the Holstein model to investigate single- 
particle excitations and phonon softening (23). However, 
it turns out that often the interaction term in Tix is cru- 
cial for a proper calculation of the required expectation 
values. Thus, usually a more involved approximation has 
been used that neglects the A dependence of the expecta- 
tion values but includes interaction effects by calculating 
the expectation values with respect to the full Hamil- 
tonian Ti. instead of Tix- In this case, the renormaliza- 
tion equations need to be solved in a self-consistent man- 
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ner because they depend on expectation values defined 
with respect to the full Hamiltonian 7i which are not 
known from the very beginning but can be determined 
from the fully renormalized (and diagonal) Hamiltonian 
H = limA-,0 Ux. 

There exist two ways to calculate expectation values 
of the full Hamiltonian from the renormalized Hamilto- 
nian. The first one is based on the free energy that can 
be calculated either from the original model Ti or the 
renormalized Hamiltonian Ti, 



F 



InTre 



-pH 



In Tr e 



because Ti. is obtained from T-L by unitary transforma- 
tions. The desired expectation values can then be de- 
termined from the free energy by functional derivatives. 
This approach has advantages as long as the derivatives 
can be evaluated analytically as, for example, in Refs.[20l 
andHH. 

The second way to calculate expectation values of the 
full Hamiltonian employs unitarity for any operator vari- 
able A, 



{A) = 



Tr(^, 



Tr 



(At 



Tre-/5« 



Tr e-P'^ 



where we defined A = lim^^o-^A- Thus, additional 
renormalization equations need to be derived for the re- 
quired operator variables A\ where the same sequence of 
unitary transformations has to be applied to the operator 
variable A as to the Hamiltonian Ti. 



where aj^^ and are given by linear combinations of 



the original fermionic operators cj^ and 



Pi 



km 



/Li + «k 4„ 



km 



km ' 



(2.23) 
(2.24) 



2 V W^k 
2 V l^k 



^^k 



Here, we defined M^k = ^^(ek-e/) +4|ykp, and the 
eigenvalues of Ti are given by 

^t'' = ^ ± \W.. (2.25) 



In the following, we want to apply the PRM as intro- 
duced above to the Fano-Anderson model (|2.2ip where 
we mainly use the formulation of Ref. ((2ll). The goal is 
to integrate out the hybridization term Tii so that we fi- 
nally obtain an effectively free model. Therefore, having 
in mind the exact solution of the model, we make the 
following renormalization ansatz: 



Ti-o.x 



no,x+ni,x, (2.26) 

(^k,A /km/km + ^k,X ^Ln'^km^ , 



k,m 



i/km^ J 



E. Example: Fano-Anderson model 

In this subsection we want to illustrate the PRM dis- 
cussed above by considering an exactly solvable model, 
namely the Fano-Anderson model (2 25l ). 



Ti. — Ti-Q + Til , 

k,m 

= E^k(/Lck. 

k,m 



(2.21) 



"km/' 



km 



The Hamiltonian (|2.2ip describes dispersion-less / elec- 
trons interacting with conduction electrons where all cor- 
relation effects are neglected, k denotes the wave vector, 
and the one-particle energies are measured with respect 
to the chemical potential. Both types of electrons are 
assumed to have the same orbital index m with values 
1, . . . jiyj. The model (|2.2ip is easily diagonalized. 



k.m 



,"kT, 



k.m 



/3LAm, (2.22) 



Note that Vk.A includes a cutoff function in order to en- 
sure that the requirement QxTi-x = is fulfilled. 

In the next step we want to eliminate excitations with 
energies within the energy shell between A and A — AA 
by means of an unitary transformation similar to (j2.17p . 
By inspecting the perturbation expansion corresponding 
to subsection lll.Bl the generator of the unitary transfor- 
mation must have the following form: 



X 



A,AA 



^Ak(A,AA) (/t„ck„ 



"k„j/km^ , 



(2.27) 



where the parameters Ak{X, AA) need to be properly de- 
termined so that Eq. (|2.18p is fulfilled. To evaluate the 
transformation (|2.17p . we now consider the transforma- 
tions of the operators appearing in the renormalization 
ansatz (|2.26p . For example, we obtain 



"-km^km L 



'^km'-km — 

= i {cos [2Ak(A, AA)] - 1} f cLckm - /L/km 

"km/km ) • 



sin[2Ak(A,AA)] (/^„Ck„ 
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Here it is important to notice that due to the fermionic 
anti-commutator relations the different k are not coupled 
with each other. Very similar transformations can also be 

found for fl^fkm and (/k„Ckm + 4m /km) ■ Inserting 
these transformations into (|2.17p leads to the following 
renormalization equations: 



= {cos [2Ak(A, AA)] - 1} [el^ ~ eQ 
+ Vu,A sin [2Au(A, AA)] , 



quantities like quasiparticle energies [compare (|2.25p with 
Eqs. (|2.32p and (|2.33p ] or expectation values. However, 
there is also an important difference between the two ap- 
proaches: Whereas the eigenmodes aj^.^ and /3j^,^j of the 
diagonalized Hamiltonian (|2.22p change there character 
as function of the wave vector k [compare p.23p and 
(I2.24p ]. the operators /^^ and of H remain /-like 
and c-like for all k values. In return, the quasi-particle 
energies e^^ and show a steplike behavior as function of 
k at £/ — Ek = so that the deviations from the original 
one-particle energies £/ and ek remain relatively small 
for all k values. 



^k,(A-AA) - £k,A - - {fixx-AX)~^L,\)- (2-29) 

Now we need to determine the parameters j4k(A, AA). 
For this purpose we employ the condition (|2.18p : First, 
from Qa^a = wc conclude Vk,A = 0k,AVk, where we 

have defined 6kA = Q (-^ ~ kk a ~ aI) • Moreover, 
from Q(A-AA)^(A-AA) = we find 

tan[2Ak(A, AA)] = (2.30) 

r 1 A 

= [1 - 0k(A-AA)J ©kA J _' ^ 

^k,A ^k,A 

which shows that also ^k(A, AA) contains the cutoff 
factor Ok, A- Note that in the expression (|2.30p the 
low excitation-energy part of the generator was chosen 
to be zero P(a-aa)-''^a,aa = 0. As one can see from 
Eqs. (|2.28p - (|2.30p . the renormalization of the parame- 
ters of a given k is not affected by other k values. Fur- 
thermore, it is important to notice that j^k a ~ ^IcaI — 

l^k (A-AA) ~ ^ScA-aaI- Consequently, each k value is 
renormalized only once during the renormalization proce- 
dure eliminating excitations from large to small A values. 
Such a steplike renormalization allows an easy solution 
of the renormalization equations (|2.28p - (|2.30p where A is 
replaced by the cutoff A of the original model and we 
set A — AA = 0. Here, one needs to consider that the 
parameter v4k changes its sign if the difference £/ — Ek 
changes its sign. Thus, we find the following renormal- 
ized Hamiltonian 

^ 1"?0'^^ " H (^k/km/km + SkcLckm) , 

k,7n 

(2.31) 

where the renormalized energies are given by 

sL = ^ + ^i^fc^I^k, (2.32) 

el ^ '-l±^ - ^g"(^;~^-V k. (2.33) 

The results of the renormalization and the diagonal- 
ization are completely comparable for physical accessible 



F. Generalized generator of the unitary transformation 

As already mentioned in subsection III.Di the low- 
energetic excitations included in the generator Xa.aa of 
the unitary transformation (|2.17p can be chosen arbitrar- 
ily, i.e. P(A-AA)-^A AA is not determined by the condition 

(EUl). 

In the previous subsection an approach of "minimal" 
transformations has been applied to the Fano-Anderson 
model where P(a-aa) A^a,aa is set to zero. However, in 
the following we want to demonstrate that it is also pos- 
sible to take advantage of this freedom to choose the gen- 
erator Xa.aa and to derive a continuous version of the 
PRM. As it will turn out in Sec. IIII.BI thc PRM can also 
be connected to Wegner's flow equation method H). 

By allowing a nonzero part P(a-aa) ArA,AA 7^ the 
generator Xa,aa of the unitary transformation (|2.17p can 
be written as follows 

ATa.AA = P(A-AA)-'^A,AA + Q(A-AA)-^A,AA (2.34) 

Here the part Q(a-aa)-'^a,aa ensures that Eq. (|2.18p . 
Q(A-AA)'^(A-AA) = 0; is fulfilled. Note however, one 
may also choose the remaining part P(A-AA)ArA,AA in 
such a way that it almost completely integrates out all 
the interactions before the cutoff energy A approaches 
their corresponding transition energies. 

As it will be discussed in Sec. IIIII in more detail, the 
flow equation method (Q) and the PRM (in its minimal 
form) take advantage of the freedom to chose the genera- 
tor of the unitary transformation in a very different way. 
In the PRM, the low transition-energy projection part of 
the generator, PaATa, is set to zero for convenience. The 
flow equation approach instead uses exactly this part to 
eliminate the interaction. 

Even though the PRM resembles the similarity trans- 
formation (6; 7) and Wegner's flow equation method (d) 
in some aspects there is an important difference: The lat- 
ter two methods start from continuous transformations 
in differential form. This has the advantage that one can 
use available computer subroutines to solve the differen- 
tial flow equations. In contrast, the PRM is based on 
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discrete transformations which lead to coupled difference 
equations. The advantage of the PRM is to provides a di- 
rect link to perturbation theory (as already discussed in 
subsection III . B I) . Moreover, the stepwise renormalization 
of the PRM allows a unified treatment on both sides of a 
quantum phase transition (see for example Sec. lVI|) which 
seems not to be possible in the flow equation method. 
However, as we show in the following the idea of continu- 
ous unitary transformations can also be implemented in 
the framework of the PRM. 



G. Fano-Anderson model revisited 

Now we want to demonstrate that the freedom in 
choosing the generator of the unitary transformation can 
be employed in order to derive a continuous renormal- 
ization scheme within the framework of the PRM. As an 
example we again discuss the Fano-Anderson model. 

As already discussed, the part P(a-aa)-'^a.aa of the 
generator AT^.aa of the unitary transformation is not 
fixed by the PRM. In the former treatment of the 
Fano-Anderson model in subsection III. El we had chosen 
P(A-AA)-^A,AA — for simplicity. In the following we 
want to take advantage of this freedom in a different way. 

According to Eq. (|2.27p . the generator of the Fano- 
Anderson model is given by 



X 



\,A\ 



k,m 



A, 



:(A,AA) (/, 



-) 



where the most general form of j4k(A, AA) can be written 
as 

Ak(A,AA) = A^A, AA) Gk.A [1 - ek,A-AA] 

-fA;^(A,AA)ek,Aek.A-AA. (2.35) 

Here, the renormalization contributions related with 
P(A-AA)-^A.AA and Q(A-AA)-'«^A.AA are described by the 
parameters A^{X,A\) and A'^{X, AX), respectively. 



A possible choice for ^k(A, A A) is 



<(A,AA) 



^k,A - ek,A ) ^k,A 



^ ^ ^k,A ^k.A 



■AA. (2.36) 



Of course, there is no derivation for Eq. (|2.36p but it will 
turn out that this is indeed a reasonable choice. In par- 
ticular we will show that in the limit of small AA a rapid 
decay for the hybridization Vk^A is obtained in this way. 
Thus, the part ^k(A, AA) of the generator is not impor- 
tant anymore for the renormalization procedure and can 
be neglected in the following. In Eq. (|2.36p . k denotes 
an energy constant to ensure a dimensionless A^{X, AA). 
Note that A'^{X, AA) is chosen proportional to AA to re- 
duce the impact of the actual value of AA on the final 
results of the renormalization. 



In order to derive continuous renormalization equa- 
tions note that the parameter Ak(A, AA) is approxi- 
mately proportional to AA. By neglecting the part 
A\^{X, AX) of the generator one can rewrite Eqs. (|2.28p 
and (|2?29)) in the limit AA ^ 



de 



k,A 



dX 
dX 



= -2Vk,xak{X) 
- +2T4,Aak(A) 



(2.37) 
(2.38) 



where higher order terms have been neglected. Further- 
more, we defined 



Q!k(A) 



lim 

AA^O 



A'^{X,AX) 



AX 



(2.39) 



•k,A ^k,A 



) ^k,A 



2 • 



A similar equation can also be derived for Vk,A 



dVi 



k,A 



O^k, 



A- 



dX ^ •'^^ "•'^^ 

To solve these equations we rewrite (|2.40p 

1 dVux 



"k.A 



e! ^ - X dX 

^k,A "'k,A 



(2.40) 



(2.41) 



-/ 



and insert into (|2.39p . Using ^ ^k a 
obtain 







d 

dX 



-(£k + 4)ek,A + Vk 



'a} - (2.42) 



Eq. (|2.42p is easily integrated and leads to a quadratic 
equation for = limA^o £k a which corresponds to the 
former result (|2.33p . Moreover, ej^ is found from ^k a + 
e( e£. According to f^AO^ and ([OOll the A- 



-k,A 



dependence of Vk.A is governed by 



dlnVj,,. 
dX 



is 



k,A ^k,A^ 



■e(A-|<A-4,Al 



(2.43) 



As one can easily see from Eq. ()2.43p . 



(i) the interaction Vk^A is always renormalized to 
smaller values when the cutoff energy A is lowered. 



(ii) and at A 



-k,A ^k,A 



the renormalized cou- 



pling Vk.A vanishes, i.e. it has completely inte- 
grated out by the present choice of the generator 
Pa-aa-'^^a-aa- 
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III. RENORMALIZATION OF THE 
ELECTRON-PHONON INTERACTION 



The classical BCS-theory (1261) is essentially based on 
attractive electron-electron interactions (f27l ). It is well- 
known that such an interaction can be mediated via 
phonons coupled to the electronic system (|28l). In this 
section we want to revisit this problem because it has 
been studied (|20l : [29I : Isol ) by Wegner's flow equation 
method ([l)j by a similarity transformation pro posed by 
Glazek and Wilson (0; 0), and by the PRM JT^). There- 
fore, the electron-phonon interaction is a perfectly suited 
test case to discuss differences and similarities of the 
three methods. In this section we consider the follow- 
ing Hamiltonian 



H = XI 4^Ck^ + ^ 6]j6q (3.1) 
k,cr q 

+ H f q ^qcLc(k+q)^ + ^qC(k+q)^Cko 



k,q,a" 



which describes electrons cj^ ^ and phonons foj^ that inter- 
act with each other. 



In the following we apply a slightly modified version 
of the PRM to the electron-phonon problem (|3.ip in or- 
der to derive an effective electron-electron interaction. 
It turns out that Frohlich's transformation (28) is re- 
examined in this way. 

In IIII.BI the approach is modified in the spirit of the 
ideas developed in subsections III.FI and III.GI Thus, al- 
lowing a more continuous renormalization of the electron- 
phonon interaction we derive the result of Ref. 29 ob- 
tained by the flow equation method. 

In subsection IIII.CI a much more sophisticated scheme 
is introduced by adding a symmetry breaking field to the 
Hamiltonian so that a gap equation can be derived. The 
effective electron-electron interaction is then obtained 
by comparing with the famous BCS-gap equation. The 
strategy to introduce symmetry breaking fields turns out 
to be of general importance for the investigation of phase 
transitions within the PRM. 

Finally, the different results for the electron-phonon 
interaction p.ip are discussed in subsection IIII.D I 



A. Frohlich's transformation 

In this subsection we want to apply the PRM to the 
electron-phonon problem (|3.1[) in order to derive an ef- 
fective electron-electron interaction. Here, we start from 
the renormalization ansatz, 

H\ = Ha + Hi, A, (3.2) 

k,cr q 



'Hi, A 
.T/cl,ph 



/T/Cl,Cl 

^1,A 



/wel.ph . /T_/el,el 



J2 [fk.q.A^Lq 

+ .9k+q,-q,A bq 



k,q,a" 



-(k+q)o- 



^k,k',q,A C(k+q)CT4k'-q)<T''^''''^"^'"^' 

k,a",k' ,0"' ,q 



that was also used in Ref. |23 where the flow equation 
method was applied to the same system. Note that the 
parameters of 7ii,A contain a cutoff function in order to 
ensure that only transitions with energies smaller than A 
are included. The parameters of Tix depend on the en- 
ergy cutoff A because all transitions with energies larger 
than A have already been integrated out. However, we 
shall restrict ourselves to the second order renormaliza- 
tion contributions to Hi, a- Therefore, Ho is assumed to 
be A independent. 

In the following we want to integrate out all transi- 
tions which create or annihilate phonons, however keep- 
ing all electronic transitions. Therefore, the present cal- 
culation differs from the previous ones where all parts of 
the 'unperturbed Hamiltonian' Ho, a were subject to the 
renormalization procedure. As it turns out, the electron- 
phonon coupling will be replaced by an effective elec- 
tron-electron interaction. However, the flnal Hamilto- 
nian containing the electron-electron interaction is not 
diagonal any more as required for the standard PRM. 
Instead, we want to derive a block-diagonal Hamiltonian 
so that the renormalization approach has to be modified. 
For this purpose, we define projection operators P^'^ and 
Qf that are defined with respect to the phonon part of 
the unperturbed Hamiltonian Ho- These new projectors 
now replace those of the full unperturbed Hamiltonian. 

Thus, from Q^'^Hi,a = we conclude 5k, q, a = 
0q,A5k,q,A, where we have defined Oq,A = 6(A — Wq). 
Moreover, following Ref. [2^, the generated electron- 
electron interaction H°'^' is not considered in determin- 
ing the generator of the unitary transformation (|2.17p . 
Thus, the generator can be written as 



AA 



(3.3) 



A,q(^> [&q4aC(k+q)^ ^ ^qC(k+q)<TCk<T 

k,q,cr 



where the parameter ylk,q(A, AA) needs to be properly 
determined in the following: Corresponding to (|2.18p . 



Q 



ph 



(A-AA)^(A-AA) 







(3.4) 



must be fulfilled. 

As already discussed, the part P(a'-aa)"'^A-^-'^ '-'^ 
generator p.3p of the unitary transformation is not fixed 
by the PRM. Thus, the parameters ^k.q(A, AA) have the 
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following general form 

Ak,q(A,AA) - A' (A,AA)eqa[l-eq,A-AA] 



<,q(A,AA) eq,Aeq,A_ 



AA- 



(3.5) 

Note that both parts of ^k.ql-^, AA) include the fac- 
tor 6q,A- However, in the following P(a-aa)-^a,aa and 
v4J^q(A, AA) are set to zero for simplicity. Note that a 
different choices for A'^ ^(A, AA) will be used in the sub- 
sequent subsection. 

We restrict ourselves to second order renormalization 
contributions so that the unitary transformation (|2.17D 
can easily be evaluated where operator terms are only 
kept if they are included in the ansatz (j3.2[) . Thus, we di- 
rectly obtain difference equation for the electron-phonon 
coupling. 



5k,q,A-AA — .9k, q, A — 
= — [^k+q — £k + t^qj 

and for the effective electron-electron interaction, 



t^al Ak+q,-q(A, AA), 



^,k',q,A-AA — ^,k',q,A = 

= ^ ^k'-q,q(A, AA).gk+q,-q,A 
-^k' -q(A, AA)5k,q,A 
1 



(3.6) 



(3.7) 



(£k+q - Ek - i^q) ^k'-q,q(A, AA) Ak,q(A, AA) 



-I- - (Ek+q - Ek + UJq) Ak+q,-q(A, AA) Ak',-q(A, AA). 

Because we have set P(a-aa)-''^a,aa = 0, renormalization 
contributions only appear if the phonon energy Uq is in 
the energy shell between (A — AA) and A. Consequently, 
we find a step- like renormalization of the electron-phonon 
coupling .gk.q.A and the generated electron-electron inter- 
action Vk_k'.q.A- The parameter ylk,q(A, AA) defined in 
(|3.5p has to be chosen in such a way that ffk.q.A-AA = 
0q,A-AA5k,q,A-AA- From equation (|3.6p we obtain 



^k,q(A,AA) = 



5q 



e 



Ek — £k+q 



-UJa 



q,A 



[1-e, 



q,A-AAj 



(3.8) 



As one can see by inserting Eq. p.Sp into (|3.6p . the 
electron-phonon coupling has no /c-dependence in the 
present approximation, i.e. gk,q,\ = ffq.A. 

Now we insert Eq. (|3.8p into the renormalization equa- 
tion (|3.7p and consider the limit A ^ 0, 



Wq |5ql 



(Ck-l-q — Ek) 



14;,k'.q — lim 14;.k'.q.A 

^ A^O 



where we exactly find Frohlich's result ([2; 



(3.9) 



B. Continuous transformation 

Wegner's flow equation method (Q) was applied to the 
electron-phonon system (|3.ip in Ref. [2^ where a renor- 
malization ansatz similar to (13.21) was used. However, a 



less singular expression for the effective electron-electron 
interaction could be derived in this way. In the follow- 
ing we want to analyze how this different result can be 
understood in the framework of the PRM. 

In order to derive continuous renormalization equa- 
tions the part Pj'a-aa)"''^-'^-^^ generator of the uni- 
tary transformation is chosen to be non-zero so that now 
A'^^{X,AX) needs to be considered in Eq. (|3.5p . Fur- 
thermore, Aj^q(A, AA) can be neglected if A'^^{\,AX) 
leads to a rapid decay of the interaction terms. Thus, 
neglecting A'^ ^(A, AA) and employing the limit AA 
we obtain from Eqs. p.6p and p.7p 



_d_ 

dA' 
d 

TT^,k',q,A 

dA 



= [<Ek+q - £k + ^q] ak,q,A, (3.10) 



.9k+q,-q,A ak',-q,A 
+ 5k,q,A Q;k'+q,q,A- 



(3.11) 



Here, we introduced Q;k,q,A = liniAA-+o ^k q(A, AA)/AA. 
Again the parameter A'^ ^(A, AA) is chosen proportional 
to AA so that the third and the fourth term on the right 
side of Eq. p.7p can be neglected in the limit AA 0. 

The commonly used generator of the flow equation 
method is chosen in such a way that the matrix elements 
of the interaction, which shall be integrated out, show 
an exponential decay with respect to the flow parameter. 
Consequently, all matrix elements change continuously 
during the renormalization procedure. We adapt the idea 
of such a continuous renormalization and assume an ex- 
ponential decay for the electron-phonon interaction. 



5k,q,A 



5q exp ■ 



(Ck+q — Ek + ^^q)" 
K (A - Wq) 



e(A 



(3.12) 



where k is just a constant to ensure a dimensionless ex- 
ponent. Note that ansatz (|3.12p is inspired by the results 
of Ref. I29I. Of course, Eq. p.l2p is only useful as long 
as the considered renormalization contributions are re- 
stricted to second order in the original electron-phonon 
interaction. Note also that ansatz (|3.12p meets the basic 
requirement ([2A8| of the PRM, Q(a-aa)W(a-aa) = 0. 

Now we need to determine the parameter ak.q.A of the 
unitary transformation. For this purpose, Eq. (j3.10p is 
divided by 5k, q, a and integrated between the cutoff A > 
Wq and 00 by using Eq. (|3.12p . We find 



ak,q,A 



.9k, q, A [Ck+q — £k 
K (A - Uqf 



(3.13) 



Note that this result is equivalent to the choice for 
A'^ q(A, AA) used for the Fano- Anderson model in III.GI 
[compare with equations (|2.36p and (|2.39p ] . 

Using this solution and the ansatz (I3.12p for the 
electron-phonon coupling (7k,q,A, Eq. (13. lip is easily in- 
tegrated where the constant k is canceled. Thus, the 
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renormalized values Vk,k',q = liniA^o Vk,k',q,A can be ob- 
tained and reads 



^k,k',q 



(3.14) 



l5ql^ (£k'. 



^q) 



(ffk+q - £k' + Wq)^ + (ek'-q " ^k' " ^^q)^ 

Iffql^ (£k'-q - £k' + ^q) 

(£k+q — Ek' — ^^q) + (£k'-q " Ek' + ^^q) 

This is the final version of the effective electron-electron 
interaction after eliminating the electron-phonon interac- 
tion. Obviously, (fXTI)) differs from Frohlich's rcsuh (28) 
that had been derived above (13.91). H owever. Eq. p.l4p 
coincides with the result of Ref. |23 that had been ob- 
tained by Wegner's flow equation method (j^). 

At this point it is important to notice that the ap- 
proaches of nil. Al and IIII.Bl are based on the same renor- 
malization ansatz (|3.2p . Therefore, the different results 
are only caused by different choices for the generator. 
Due to the continuous renormalization, the electron- 
phonon coupling becomes dependent on the electronic 
one-particle energies £k so that the approach of llll.BI in- 
volves more degrees of freedom. 

The main goal of this subsection was to demonstrate 
that Wegner's flow equation method (8) can be under- 
stood within the PRM(IO), as already for the case of the 
Fano- Anderson model in the previous section. However, 
the idea of a continuous renormalization, as implemented 
here, can also be very useful for other applications. In 
this regards, the discussion line needs to be changed: 
One starts from an ansatz for the generator ATa^aa of 
the unitary transformation similar to Eqs. I|3.3p . (I3.13p . 
and demonstrates afterwards that the interaction decays 
as function of A as required. 



C. Improved renormalization scheme and BCS-gap 
equation 

So far the discussion of the electron-phonon problem 
was focused on the phonon-induced electron-electron in- 
teraction. Thus, we derived block-diagonal Hamiltonians 
with constant phonon occupation numbers within each 
block. However, in the following we want to tackle the 
electron-phonon problem (|3.ip in a different way because 
an effective phonon mediated electron-electron interac- 
tion is mainly discussed with respect to superconductiv- 
ity. The idea is to obtain the superconducting properties 
directly from the electron-phonon system. 

The goal is again to decouple the electron and the 
phonon system but now we want to derive a truly di- 
agonal renormalized Hamiltonian. For this purpose the 
PRM shall be applied to the electron-phonon system 
( 3.ip in conjunction with a Bogoliubov transformation 
( 31|) as it was done in Ref. |20|. 

Whereas the Hamiltonian p.ip is gauge invariant, a 
BCS-like Hamiltonian breaks this symmetry |26[). There- 
fore, in order to describe superconducting properties, the 



renormalized Hamiltonian should contain a symmetry 
breaking field as well so that the renormalization ansatz 
reads 



'Ho,\ + Hi, A, 



(3.15) 



J2 (^k.A 4^cLkx + K,X C-kiCkt) + Cx 



Hi,. 



k,q,cj 



'^ko-'^(k+q)o-^q + C(k+q)cr'^kcr&q 



Here, the 'fields' Ak,A and AJ^ ^ break the gauge invari- 
ance and can be interpreted as the superconducting gap 
function. The initial values for Ak,A and the energy shift 
Cx are given by those of the original model, Ak,A ~ 0, 
Ca — 0. Note that in the following the projectors Pa 
and Qa are defined as usual with respect to Tio.A and 
not only to the phonon part. Furthermore, renormaliza- 
tion contributions to electronic and phononic one-particle 
energies and to the electron-phonon coupling will be ne- 
glected for simplicity. 

At this point it is important to realize that the intro- 
duction of symmetry breaking fields is a general concept 
to study phase transitions within the PRM. The same 
approach has also been successfully applied to the Hol- 
stein model and its quantum phase transition ((3^ : [sl); 
this model will be discussed in Sec. IVII 

To perform our renormalization scheme as introduced 
in section [TT] we need to solve the eigenvalue problem of 
Ti-o.x- For this purpose we utilize the well-known Bogoli- 
ubov transformation ((31i ) and introduce new A dependent 
fermionic operators. 



4a 



* t * 

k.A^kt "k.A^-ki, 



(3.16) 



Plx = «: 



k,AC_ki +"k,ACkT, 



where the coefficients read 




(3.17) 

/4+iAk,Ar^ ' 

Hence, 7io,A can be rewritten in diagonal form, 

no,x = X ^k,A (a^akA + /3L/3kA) (3.18) 

k 

+ X (ek - £^k,A) + X t^q b^b^ + Cx 

k q 

where the fermionic excitation energies are given by 
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In the following, we restrict ourselves to second order 
renormalization contributions so that the first order of 
the generator Xa,aa of the unitary transformation is suf- 
ficient [see Eq. (|2.12p and the discussion in llLBj . Thus, 
-'^^A.AA can be written as (|2.17p . 

^A,AA = (3.19) 

= J2 ^k,q(A, AA) [ fo]j4^C(k+q)^ - ^c[i^+q)<xCk<T 
k,q,(7 



where 

^k,q(A,AA) 



ek,q(A,AA), 



(3.20) 



ek,q(A,AA) = [l-e(A-AA-|ek-ek+q 

X6 (A - |£k - £k+q +t^q|) . 



Note that the generator Xa,aa as defined in Eq. (|3.19p al- 
most completely agrees with the one used to re-examine 
Frohlich's transformation in subsection IIII.AI [see Eqs. 
(13. 3p and p.Sp ]. However, now the Q functions do not 
only refer to the phonon energies Wq but also to the elec- 
tronic one-particle energies £k because of the different 
definitions of the Pa projection operators. 

To perform the renormalization step reducing the cut- 
off from A to A — AA, one would need to express the elec- 
tronic creation and annihilation operators by the quasi- 
particle operators p.l6p . After considering the renormal- 
ization contributions, the quasi-particle operators have 
to be transformed back to the original electron opera- 
tors. However, this involved procedure is only necessary 
if we are interested in renormalization contributions be- 
yond second order perturbation theory. Therefore, here 
the symmetry breaking fields Ak,A and ^ are only gen- 
erated by the renormalization scheme but not considered 
in the evaluation of energy denominators or projection 
operators. 

Taking into account all simplifications related with sec- 
ond order perturbation theory, the unitary transforma- 
tion (|2.17p is easily evaluated where generated opera- 
tor terms are only kept if their mean-field approxima- 
tions renormalize the symmetry breaking fields, Ak,A and 
A^ A' '^^ energy shift, C\. Thus, for sufficiently small 
steps AA we obtain the following renormalization equa- 
tions 



Ak,A-AA — Ak,A = 
= 2 ^ e [A - jsk - £(k+ci) I + ^^q] 

q 

X {l - e [A - AA - |£k - £(k+q) 

|.9q|^0 [t^q- l^k - e(k+q)'"' 



(A-AA) 



^(k+q) 



E 

k 



(3.21) 

-J} 

(c-(k+q),iC(k-|-q),T) , 



''k,T'--ka 



[Aw.A- 



k.A-AA 



A' 



k.A 



(3.22) 



By summing up all difference equations between the cut- 



off A of the original model and the lower cutoff A 
one easily finds 



0, 



Ak = A 



k,A 



l5qre 



^(k+q) 



q 



X (c 

E 



|£k - £(k+q)| 
(k+q),iC(k+q),T) , 



^k,T'--k,| 



Ak- A 



k.A 



(3.23) 



(3.24) 



Here we defined Ak = liniA^o Ak,A, C = limA^o C\- 

The final Hamiltonian Ti. = limA^o Ti-x can easily be 
diagonalized by a Bogoliubov transformation and reads 
according (|3.18p 



n = 



^ ^k (5|,5k + PlPk 



(3.25) 



^(ek-^k) +5I^q^q^q 



C 



where £'k = limA^o -E^k.A, Sk = limA^oak,A, and /3k = 
limA— ►0 /3k, A. Its parameters depend on the original sys- 
tem (|3.1I) . on the initial conditions, Ak,A ~ 0, Ca = 0, 

and on expectation values ^cj^jC^j^i^ that need to be 

determined self-consistently. Following the approach of 
Ref. |20|, we consider the free energy which can be calcu- 
lated either from Ti. or from the renormalized Hamilto- 
nian Ti.. Thus, the required expectation values are easily 

4 J \-_dF_ 
Ck,TC-ka/- 

so that Eq. (|3.23p can be rewritten as 



found by functional derivatives. 



Ak = E 



2|5q|2e[a;q- 


£k — £(k-|-q) 


] 




Ek - £(k+q) 





(3.26) 



^k+q 



1 - 2/(^k^ 



Ak 



where the initial condition Ak,A — has been used. 
Eq. p.26p has the form of the famous BCS-gap equa- 
tion so that the term inside the braces {• ■ • } can be in- 
terpreted as parameter of the effective phonon induced 
electron-electron interaction, 



Vk 



-k,q 



|.gq|^B [Wq- |£k-e(k+q)|] 



|£k - £(k+q)| 



(3.27) 



which is responsible for the formation of Cooper pairs. 
Even though we have here derived an effective electron- 
electron interaction as well there is a significant difference 
to the approaches of IIII.AI and IIII.Bl In the present for- 
malism both the attractive electron-electron interaction 
and the superconducting gap function were derived in 
one step by applying the PRM to the electron-phonon 
system (j3.ip with additional symmetry breaking fields. 
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D. Discussion 

In the following we want to discuss the different 
approaches to the phonon-induced effective electron- 
electron interaction in more detail. At first we summa- 
rize the results derived above where we focus on the in- 
teraction between electrons of a Cooper pair. Frohlich's 
classical result [see Ref. i28| and Eq. (|3.9p ] reads 



Fr6hl] 



k,q 



(3.28) 



However, there is an important problem related with 
Eq. p.28p : It diverges at |ek — £(k+q)| = t^q- Thus, 
a cutoff function is introduced by hand in the classical 
BCS-theory to suppress repulsive contributions to the ef- 
fective electron-electron interaction. 

In contrast to the Frohlich interaction ()3.28|) . the re- 
sults obtained by Wegner's flow equation method (l29h . 
by similarity transformation ((sol). and by the PRM (|2Clt ) 
are less singular. 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



(e, - 8, ) / CO 



FIG. 1 Comparison of the effective electron-electron interac- 
tion obtained by the PRM (full line) and by Wegner's flow 
equations (dashed line). Here, the dimensionless quantity 
iikq = — ^ has been introduced. 



T rLonz/Wogn 

k,-k,q 



T/Miclkc 

''k.-k.q.A 



T T-Hiibsch/Bcckc 

''^k,-k,q,A 



ISql^^q 



(Sk+q-fk) 
|5ql^0(kk+q-ek| 



(3.29) 
A) 



|£k+q - Ekl 



|.gq| Q (^q - kk+q - £k | 
kk+q - £k| + ^q 



(3.30) 



(3.31) 



(Note that Eqs. (|3.29p and (|3.3ip have already been de- 
rived above, compare with (|3.14p and (|3.27p . The A 
dependence of the electronic and phononic one-particle 
energies are suppressed in p.30|l for simplicity.) All 
three results for the effective phonon-mediated electron- 
electron interaction are never repulsive as long as > 
is fulfilled. 



At first we want to discuss Mielke's result fSCr), an effec- 
tive electron-electron interaction p.30p that depends on 
the energy cutoff A. As Wegner's flow equation method 
(H), the used similarity transformation H; 0) is based 
on continuous unitary transformations and leads to dif- 
ferential equations for the parameters of the Hamilto- 
nian. However, like the PRM, the similarity transfor- 
mation leads to a band-diagonal structure of the renor- 
malized Hamiltonian with respect to the eigenenergies of 
the unperturbed Hamiltonian whereas the flow equation 
method generates block-diagonal Hamiltonians. 

Mielke derived the phonon-mediated electron-electron 
interaction p.30p by eliminating excitations with ener- 
gies larger than A where excitation energies are measured 
with respect to the unperturbed Hamiltonian consisting 
of both electronic and bosonic degrees of freedom. The 
obtained effective interaction becomes A independent for 
the Einstein model (of dispersion-less phonons) if A is 



chosen smaller than the phonon frequency uq- For this 
case Mielke's result p.30p is very similar to ours ()3.3ip 
obtained by the PRM with symmetry-breaking fields. 
However, in contrast to our result (|3.3ip . the cutoff func- 
tion 9 (wq — |ek-i-q — Ekl) is absent in (|3.30p . This dif- 
ference might be related with different choices for the 
generator of the unitary transformation in the two meth- 
ods but could also be caused by a systematic problem in 
Mielke's approach: Setting A = 0, the final renormalized 
Hamiltonian contains non-diagonal terms with respect to 
the used unperturbed Hamiltonian. This seems to con- 
tradict a basic premise of the similarity transformation. 



Lenz and Wegner (|29[ ) applied the fiow equation 
method to the electron-phonon problem as discussed here 
and obtained an effective electron-electron interaction as 
shown in Eq. ()3.29p . As one can see in Fig.[l] their result 
is quite similar to ours p.3ip derived using the PRM as 
long as LOq > l^k+q — £k| is fulfilled. However, in contrast 
to our result p.3ip . the interaction (|3.29p remains finite 



even for 



< |ek4 



ek|- Probably, this difference is 



caused by the different choices for the generator of the 
unitary transformation that also require different approx- 
imations in order to obtain closed sets of renormalization 
equations. 



IV. HEAVY-FERIVIION BEHAVIOR IN THE PERIODIC 
ANDERSON MODEL 

The periodic Anderson model (PAM) is considered to 
be the basic microscopic model for the theoretical inves- 
tigation of heavy- fermion (HF) systems IH). It describes 
localized, strongly correlated / electrons interacting with 
itinerant conduction electrons. Here we focus on the limit 
of infinitely large Coulomb repulsion on / sites so that 
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the Hamiltonian of the PAM can be written as 

n = Ho+Hi, (4.1) 



The one-particle energies £/ and Ck, and, as a simph- 
fication, both types of electrons have the same angular 
momentum index m — 1 . . .Vf. The Hubbard operators, 

ftn ^ fim ~ fimfi-rfi-)' 

take into account the infinitely large local Coulomb re- 
pulsion and only allow either empty or singly occupied / 
sites. 

The PRM has already been applied to the PAM in 
Ref. [21I : [2^ where approximations have been employed 
that allow to map the renormalization equations of the 
PAM onto those of the uncorrelated Fano- Anderson 
model (see subsection III.Ep . Thus, HF behavior and a 
possible valence transition between mixed and integral 
valent states could be studied. However, the approach 
of Refs. lU has a significant disadvantage: the renor- 
malization of the one-particle energies show as function 
the cutoff A a steplike behavior that leads to serious prob- 
lems in the (numerical) evaluation. Therefore, a constant 
renormalized / energy had to be chosen for all values of 
the energy cutoff A to ensure a continuous behavior of 
the one-particle energies as required for physical reasons. 

In the following we modify the approach of Refs.[2ll:[23 
to ensure a more continuous renormalization of all pa- 
rameters of the Hamiltonian. For this purpose, the ideas 
of III.GI and IIII.Bl are transferred to the PAA4. However, 
to explore all features of this continuous approach is be- 
yond the scope of this review, we re-derive the analytical 
solution of Ref. instead. 



A. Renormalization ansatz 

Much of the physics of the PAM (|4.ip can be under- 
stood in terms of an effective uncorrelated model that 
consists of two non-interacting fermionic quasi-particle 
bands. Various theoretical approaches have been used to 
generate such effective Hamiltonians; the most popular 
among them is the slave-boson mean- field (SB) theory 
(fill : [12 ). However, as discussed in Ref. [22, such ap- 
proaches do not prevent from unphysical multiple occu- 
pation of / sites and are therefore restricted to heavy- 
fermion like solutions. [The SB solutions break down if 
the original / level £/ is located too far below the Fermi 
level or if the hybridization between / and conduction 
electrons becomes too weak 1351) .1 



To reliably prevent the system from unphysical states 
with multiple occupations of / sites we here follow Ref. [13 
and start from a renormalization ansatz that keeps the 
Hubbard operators during the whole renormalization 
procedure. 



Hoa+Wia, (4.2) 



k.m 



k, m 



NL 



k.m 

k,m 



Eq. (|4.2p is obtained after all excitations between eigen- 
states of 7io,A with transition energies larger than the 
cutoff A have been eliminated, i.e. QxTix = holds. 
Furthermore, we introduced Fourier transformed Hub- 
bard operators, 



•> km 



1 



ik-Ri 



The A dependencies of the parameters are caused by the 
renormalization procedure. Note that Vk.A includes a 
cutoff function in order to ensure that the requirement 
Qa'^a = is fulfilled. Furthermore, an additional energy 
shift E\ and direct hopping between / sites. 



ik(Ri-Rj) 



have been generated. Finally, we need the initial param- 
eter values of the original model (with cutoff A) to fully 
determine the renormalization, 

e/,A = Sf, Ak,A = 0, ek,A = ^k, Ea = 0, (4.3) 

Vk,A = Vk. 

To implement our PRM scheme we also need the com- 
mutator of the unperturbed part 7io,A of the A depen- 
dent Hamiltonian Tix with the interaction Tii^x (in the 
present case the hybridization between / and conduction 
electrons). To shorten the notation we here introduce the 
(unperturbed) Liouville operator Lo,a that is defined as 
Lo,A-4 — [Ho,x,-A] for any operator A. Because of the 
correlations included in the Hubbard operators fl^^, the 
required commutator relation can not be calculated ex- 
actly and additional approximations are necessary. Here, 
the one-particle operators and cj^^ are considered as 
approximative eigenoperators of Lq^a so that we obtain 



I^O.A /k„Ck„ 



(£/,A + £'Aka 



^k,A)/k,i 



Ckri 



Here we introduced the local / energy, 
£/,A = e/,A - -DAa, 



(4.4) 



(4.5) 
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the averaged / dispersion, Aa = jj^k^k,X: and de- 
fined D = 1 - (nf) + {n{)/vf. Note that the factors D 
in Eqs. (|4.4p and (|4.5p are caused by the Hubbard oper- 
ators /j^^ where a factorization approximation has been 
employed. 

To ensure that QaT^a = is fulfilled by g^]), the 
hybridization matrix elements must include an additional 
6-function, Vk_A — 9(k, A)Vk,A, where we have defined 

e(k,A) = e(A- |e/,A + Z?Ak,A-ek,A|). 

B. Generator of the unitary transformation 

In order to derive the renormalization equations for the 
parameters of Tix we have to consider the unitary trans- 
formation to eliminate excitations within the energy shell 
between A — A A and A. Corresponding to Eq. (|2.17p . such 
a unitary transformation is determined by its generator 
^A,AA- As in Ref . we use an ansatz that is motivated 
by perturbation theory [see Eq. (|2.12p ]. 

Xa,aa - 5]Ak(A,AA) (/Lck™ - . 

k.m 

(4.6) 

The parameter ^k(A, AA) of the generator Xa^aa 
needs to be chosen in such a way that Eq. (|2.18p . 
Q(A-AA)W(A-AA) — 0, is fulfilled. However, as already 
discussed before, this requirement only determines the 
part Q(A-AA)-'^A,AA of the generator (|4.6p of the unitary 
transformation whereas P(a-aa)^a,aa can be chosen ar- 
bitrarily. Thus, P(a-aa)-'^a,aa = is usually chosen 
to perform the minimal transformation to match the re- 
quirement (|2.18p . In this way, the impact of approxi- 
mations necessary for every renormalization step can be 
minimized. 

On the other hand, the approach of "minimal" ap- 
proximations can also lead to some problems if a step- 
like renormalization behavior for the parameter of the 
Hamiltonian is found. This is the case for the PRM ap- 
proach of Refs. [21]; [22 where a constant renormalized / 
energy e / have been used for all cutoff values A to ensure 
a continuous behavior of the one-particle energies as re- 
quired for physical reasons. Therefore, in the following 
P(A-AA)^A,AA shall again be chosen non-zero in order to 
ensure a more continuous renormalization of all param- 
eters of the Hamiltonian. In close analogy to subsection 
III.Gi we choose a proper generator ^j^(A, AA) ~ AA, 
not yet specified, which almost completely integrates out 
interactions before the cutoff energy A approaches their 
corresponding transition energies. In the limit of small 
AA, we again expect an exponential decay for the hy- 
bridization Vk,A in this way. 

C. Renormalization equations 

In comparison to the approach of Refs. [U [2^ the 
derivation of the renormalization equation is simplified: 



Having in mind j4J^(A, AA) ^ AA, where AA is a small 
quantity, we can restrict ourselves to first order renormal- 
ization contributions and neglect the A'^{\,AX) part of 
^A.AA altogether. Thus, eliminating excitations within 
the energy shell between A — A A and A, the renormal- 
ized Hamiltonian 7i(A-AA) can be calculated based on 

Eq. mm . 

To derive the renormalization equations for the param- 
eters of the Hamiltonian, we compare the coefficients of 
the different operator terms in the renormalization ansatz 
(|4.2p at cutoff A — AA and in the explicitly evaluated 
Eq. (|2.19p . Thus, based on similar approximations as 
the approach of Refs. and '22, we obtain the following 
equations: 

ek,A-AA-ek,A = -2i?A;;(A,AA)Vk,A (4.7) 

Ak,A-AA - Ak,A = - [ek,A-AA ~ , (4.8) 

e/,A-AA-e/,A = (4.9) 

k 

X {l + (i^/-l)(cLckm)} 

+ AA) <(A, AA) (Ak,A - Aa) 

k 

X (/Lckm+h.C.^, 
Vk,A-AA - ^k,A = (4.10) 

= -A;^(A,AA) [e/,A + i?(Ak,A- Aa) -ek,A] 

£^(A-AA)-i?A = - A^(n.f ) [e/,A-AA - e/,A] (4.11) 

(^f) r 1 
jy- 2^ [£k,A-AA - ek,Aj ■ 

k 

Here, the condition Vk,A-AA = 6(k, A — AA)Vk,A-AA 
has to be fulfilled. Note that higher order terms in 
these equations have been evaluated in Refs. [2l| and 
I22 for the case that the generator Xx^ax was fixed by 

Q(A-AA)-''^A,AA- 

In deriving the renormalization equations (|4.7p - (|4.1ip 
a factorization approximation has been employed in or- 
der to trace back all terms to operators appearing in the 
renormalization ansatz (|4.2p . Thus, the renormalization 
equations still depend on expectation values which have 
to be determined simultaneously. Following the approach 
of Ref. i22, we neglect the A dependency of all expecta- 
tion values and calculate them with respect to the full 
Hamiltonian H. As discussed in subsection III.Dl there 
are two strategies to obtain such expectation values: The 
first one is based on the free energy which we will use 
later for the analytical solution in IIV.DI However, the 
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evaluation of the free energy is complicated as long as 
the renormalized Hamiltonian contains Hubbard opera- 
tors /km- Thus, here it would be more convenient to use 
the second strategy to calculate expectation values and 
to derive renormalization equations for additional oper- 
ator expressions (see Refs. [21I and [22I for more details). 
However, such involved approach is only needed in case of 
a numerical treatment of the renormalization equations 
which will be discussed below. 

The further calculations can be simplified by consid- 
ering the limit AA — > and to transform the difference 
equations (|4.7p - (|4.1ip into differential equations. For 
this purpose we define 



a.(A) = hm 



AA^O AA 



so that we obtain 

d£k,A 



dA 



= 2i?ak(A)li,A, 



dA 



k,A 



1 de 



k,A 



dA 



D dA 



(4.12) 



(4.13) 



(4.14) 



de 



/■A 



dA 



1 1 

DN 



E{i + (^/-i)(^Lck™)} 



de 



k,A 



dA 



N 



^e(k,A)ak(A) (Ak,A- A> 



dA 
dA 



X (/t,„ck™ + h.c.), (4.15) 
= [e/,A + i?(Ak,A- Aa) -ek,A] ak(A), (4.16) 



- N{h 



def.A (ni) v-^ de 



/\ae/,A 



dA Z? ^ dA 

k 



-k,A 



(4.17) 



D. Analytical solution 

In the following, we concentrate on an analytical so- 
lution of the renormalization equations (|4.13p - (|4.17p by 
assuming a A independent energy of the / electrons. 
The aim is to demonstrate that the analytical solution 
of Ref. I21I can also be derived from the renormalization 
equations (|i?T^ - (|iT7)) or likewise (|i?71) - (liTT)) obtained 
here. In particular, we want to derive an analytical so- 
lution that describes HF behavior. As in Ref. [2l|, we use 
the following approximations: 

(i) All expectation values (which appear due to the 
employed factorization approximation) are consid- 
ered as independent from the renormalization pa- 
rameter A and are calculated with respect to the 
full Hamiltonian Ti.. 

(ii) As mentioned, the A dependence of the renor- 
malized / level is neglected and we approximate 
Cf^x — DA\ if to decouple the renormalization 



of the different k values. Note that such a renormal- 
ized / energy is also used from the very beginning 
in the SB theory. 

(iii) To obtain the analytical solution of Ref. [21I we set 
jj Ak = for further simplification. 

(iv) The Hubbard operators are replaced by usual 
fermionic operators where we employ 

E'^km/km = E -^km/k™ ^^^d 



(/km/kni^ 



NL V /NL 



Thus, on a mean- field level, the system is prevented 
from generating unphysical states but a multiple 
occupation of / sites is not completely suppressed 
by this approximation. Therefore, we can only ob- 
tain useful results as long as only very few / type 
states below the Fermi level are occupied. 

It turns out that the analytical solution of Ref. [Ill is 
obtained if the approximations (i)-(iii) are applied to the 
renormalization equations (I4.13p - (l4.17p . 

Employing approximation (iv), the desired renormal- 
ized Hamiltonian Ti. = liniA^o ^a is a free system consist- 
ing of two non-interacting fermionic quasi-particle bands. 



n = E^^k4 



km'^km 



(4.18) 



k.m 



E(£/+^Ak) fl^hrn + E. 
k,m 



Eqs. (|4.14p and (|4.1ip can be easily integrated between 
A = and the cutoff A of the original model. 



A. 



D 



[ek - £kj , 



(4.19) 



/\ [f . _ . .1 ^ ^ l/f.f\ J2 [£k - £k] 



E = -N{h{)[ef^ef] + 



D 



-N{h{) [sf^Sf], 



(4.20) 



where approximation (iii) has been used. The equation 
()4.13p can also be solved if the renormalizations of the 
different k values are decoupled from each other by ap- 
proximations (i) and (ii). Thus, Eq. (|4.16p can be rewrit- 
ten as 



ak(A) = — 



1 



k,A 



e/+£k-2ek,A dA 
and inserted into (|4.13p so that we obtain 
d 







dA 



{elx - (£/ + £k) £k,A + DV^x} . (4.21) 



Eq. (|4.2ip can easily be integrated and a quadratic equa- 
tion for Ek = lim^^o £k,A is obtained. Our recent work 
on the PAM m il) has shown that the quasi-particles 
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in the final Hamiltonian (|4.18p do not change their (c 
or /) character as function of the wave vector k. There- 
fore, £k jumps between the two solutions of the obtained 
quadratic equation in order to minimize its deviations 
from the original £k, 



- £/ + £k sgn(£/ - £k) 

£k = 7; 7; ^k, 



M^k = 



- iff 



(4.22) 
(4.23) 



The second quasi-particle band is given by 

~ , nA £/ + £k , sgn(£/ ~ £k) 
Wk := £/ + DAk = 1 Wk. 

(4.24) 

Thus, we have obtained the same effective Hamiltonian 
(|4.18p and the same quasi-particle energies (|4.22p and 
(|4:24l) as found in Ref.M 

Finally, we need to determine the renormalizcd / en- 
ergy if and the expectation values. Because the fi- 
nal renormalizcd Hamiltonian (j4.18p consists of non- 
interacting fermionic quasi-particles, it is straightforward 
to calculate all desired quantities from the free energy as 
it was done in Ref. 21- Because the effective model Ti 
is connected with the original Hamiltonian Ti. by an uni- 
tary transformation the free energy can also be calculated 
from 7i, 

F = -^InTre"'^". 

P 

The expectation value of the / occupation is found from 
the free energy by functional derivative. 



(n!) 



1 OF _ I I dn\ 



(4.25) 



n 



Thus, we finally obtain a relation of the following struc- 
ture 




(4.26) 



In the cases of mixed valence and heavy Fcrmion behav- 
ior the derivatives in Eq. (|4.26p are non-zero so that both 
brace expressions can be set equal to zero to find equa- 
tions of self-consistency for the renormalizcd / level and 
the averaged / occupation number. 



N 



-+Sgn(£,-£k)^^ 



(4.27) 



k 



- sgn(£k - if 



Vf - 1 



£k - £/ 
21¥k 



^sgn(£7-£k)/(£k)^ (4.28) 



k 



+ ^V^E^g<^k-£/)/(^k) 




FIG. 2 (Color online) /-electron occupation number = 
{n{ ) as function of the bare energy e / for an one-dimensional 
lattice with 10000 sites for two cases: i) the total particle 
occupation n^ot = -\-n'^ = 1.75 is fixed (in red) and (ii) the 
chemical potential (in green) is fixed. Moreover, Uf = 2, 
V = 0.10(4f) and the temperature T = 



These equations are quite similar to the results of the 
SB theory (fT3 ). In particular, the limit j/y — > cx3 of 
Eqs. ((427|) and (|4?28)) leads to the SB equations. Note 
that expectation values (cj^^Ckm) and (/^,„Ckm, + h.c.) 
can be calculated similar to Eq. (|4.25p . see Ref. |2l| for 
details. 



E. Numerical solution 

Note that for the analytical solution in the preceed- 
ing subsection an explicit expression for the generator 
v4"(A, AA), was not needed. The reason was that a 
A independent / electron energy £/_> was assumed in 
close analogy to what is done in the well known slave 
boson mean field approach for the periodic Anderson 
model. For an improved treatment an explicit expres- 
sion for ^k(A, AA) should be used. Following the discus- 
sion in subsection III.GI we make the following ansatz for 
^i;(A,AA) 



<(A,AA) = 

(e/,A + £'(Ak,A 



£k,A) Vk,. 



K [A - |e/,A + £> (Ak,A 



(4.29) 
AA. 



£k,A 



In the limit of small AA, we again expect an expo- 
nential decay for the hybridization Vk^A in this way. In 
Eq. (|4.29|) . K denotes an energy constant to ensure a di- 
mensionless A(^(A, AA). Note that AJ^(A, AA) is chosen 
proportional to AA to reduce the impact of the actual 
value of AA on the final results of the renormalization. 
Using (|4.12|) and (|4.29[) the basic renormalization equa- 
tions (|4.13p - (|4.17p was solved numerically in Ref. |36 



FIG. [2] shows the / occupation 



i{) as func- 



tion of the bare / energy Sf at degeneracy Vf = 2 for 
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two cases, (i) for fixed total particle occupation ti^qi = 
+ n'^ — 1.75 (in red) and (ii) for fixed chemical poten- 
tial fi (in green). Here, = (l/N) I]k,<T(4,<TCk,<T) is the 
conduction electron occupation. For the first case the 
result from the PRM approach shows a rather smooth 
decay from the integer valence region with — 1, when 
£/ is located far below the Fermi level, to an empty state 
with no / electrons — 0, when £/ is far above the 
Fermi level (black line). Note that this analytical PRM 
result almost completely agrees with the result from re- 
cent DMRG calculations from Ref. 37 for the same pa- 
rameter values. For comparison, the figure also contains a 
curve obtained from the PRM approach when the chemi- 
cal potential /i instead of «tot fixed in the calculation 
(red curve). Note that in this case as function of £/ 
shows an abrupt change from an completely filled to an 
empty / state. Obviously the latter behavior can easily 
be understood as change of the / charge when e / crosses 
the fixed chemical potential. In contrast, for fixed total 
occupation n^Qf^ the Fermi level is shifted upwards, when 
the / level is partially depleted when e / comes closer to 
the Fermi level. For details we refer to ReflSGl 

V. CROSSOVER BEHAVIOR IN THE METALLIC 
ONE-DIMENSIONAL HOLSTEIN MODEL 

In this section we discuss the one-dimensional Holstein 
model. As is well known, this model shows a quantum 
phase transition between a metallic and a charge ordered 
state as function of the electron-phonon coupling. In the 
present section we restrict ourselves to the metallic state. 

Let us start with the Hamiltonian of the one- 
dimensional Holstein model of spinless fermions (HM) 
which reads, 

n = -t^(clcj + h.c.) + ^oE ^l^^ (5.1) 

i 

This model is perhaps the simplest realization of an 
electron-phonon (EP) system and describes the interac- 
tion between the local electron density Ui = c\ci and 
dispersion-less phonons with frequency loq. Here, the 
cj {b\) denote creation operators of electrons (phonons), 
and the summation (i, j) runs over all pairs of neighbor- 
ing lattice sites. With increasing EP coupling g, the HM 
undergoes the quantum-phase transition from a metallic 
to a charge-ordered insulating state. At half-filling, the 
insulating state of the HM is a dimerized Peierls phase. 

Because the HM is not exactly solvable, a number of 
different analytical and numerical methods have been 
applied: strong couphng expansions HI), Monte Carlo 
simulations Issl : Isol ) , variational (|40| ) and renormaliza- 
tion group (l4l[) approaches, exact diagonalization (ED) 
techniques (|42r). density matrix renormalization group 
dH; HI Hi) and dynamical mean-field theory (DMFT) 



(|46h . However, most of these approaches are restricted 
in their application, and the infinite phononic Hilbert 
space (even for finite systems) demands the application 
of truncation schemes in numerical methods or involved 
reduction procedures. 

The PRM represents an alternative analytical ap- 
proach. In the following the PRM is applied to the HM 
where we mainly follow Refs. 32, and 33. Here we focus 
on the investigation of the change of physical proper- 
ties by passing from the adiabatic to the anti-adiabatic 
limit. Furthermore, we discuss electronic and phononic 
quasi-particle energies as well as the impact of the system 
filling. 

A. Metallic solutions 

For the metallic phase of the HM a very simple renor- 
malization scheme is sufficient where only the electronic 
and phononic one- par ticle energies are renormalized. 

Following Refs. |23| and [H, we make the following 
ansatz for the renormalized Hamiltonian 

"Hx = Wo,A+'^i,A, (5-2) 

Wo, A = £k,\c\ck + ujq,xblbg + Ex, 

k q 
k,q 

Here, all excitations with energies larger than a given cut- 
off A arc thought to be integrated out. Moreover, we have 
defined Ofc^g, A ^ Q{X - \u!q,x + £k,x - £k+q,x\)- Note that 
Fourier-transformed one-particle operators have been 
used for convenience. Next, all transitions within the 
energy shell between A — AA and A will be removed by 
use of a unitary transformation (Eq. (|2.17p ). 

Hix-AX) = e^^-^^ Hx e-""--^^ , (5.3) 

where the following ansatz is made for the generator 
^A.AA of the transformation 

(5.4) 

The part P(a-aa)-'^a.aa has been set equal to zero. 
Therefore Ak.q(X, AX) reads 

Ak^q{X,AX) = A'kg{X,AX)ek.q^x[l-'dk.q.x-Ax]. 

As before, the ansatz (|5.4p is suggested by the form of 
the first order expression (|2.20|1 of the generator X\^ax- 
Later, the coefficients A'l^ ^(A, AA) will be fixed in a way 
that Q(A-AA)W(A-AA) =' is fulfilled, so that 7i(A-AA) 
contains no transitions larger than the new cutoff A — AA. 

By evaluating (|5.3p . terms with four fermionic and 
bosonic one-particle operators and higher order terms 
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are generated. In order to restrict the renormalization 
scheme to the terms included in the ansatz (|5.2p . a fac- 
torization approximation has to be employed, 



b%clck 



-(4cfc)(cLgCfe_,), 



In this way, it is possible to sum up the series expansion 
from transformation ()5.3p . 

The parameters A'^. ^(A, AA) as well as the renormal- 
ization equations for Ek^x, ^q,x, 9k,q,\^ and Ex can be 
found by comparing the final result obtained from the ex- 
plicit evaluation of the unitary transformation (|5.3p with 
the renormalization ansatz ()5.2|) . where A is replaced by 
A — AA. The result is given in Ref. i23i. It can be fur- 
ther simplified in the thermodynamic limit N ^ oo. By 
expanding the renormalization equations from Ref. 23 in 
powers of g, one finds that only terms of quadratic or 
linear order in g survive. The final equations read 



efe,(A-AA) — Efe.A — 

/VZ^V 9 k + q) u;q^^+£^ ,^^ek + q,X 



(5.5) 



9 



k—q 



1) 



g^efc_,,g(A,AA) 

'-^q,X + £k-q,X — £fc,A ' 



^q.iX-AX) - 
k 



(5.6) 



'■k+q) 



ff^9fc,,(A,AA) 

Wq,A + Efe.A — £fe+g,A 



"•fc = (4cfe)> (^g^g); and efc,g(A, AA) = 



where 

0fc,g,A [1 — 0A;,g,A-AAj- 

Note that the renormalization equations still depend 
on unknown expectation values (c^Cfc) and (b^qbq) which 
follow from the factorization approximation. Following 
Ref. m, they are best evaluated with respect to the full 
Hamiltonian Ti.. 

Exploiting (A) = limx^o{Ax)Hxi derive additional 
renormalization equations for the fcrmionic and bosonic 
one-particle operators, cj, and They have the follow- 
ing form according to Refs. [lH and[33l 

4,A = "fc,A 4 + E (^'=^9^^ 4+g^9 + 7fc,g.A C^g^j) , 

(5.7) 

^Ix = <l^qAbl+Vq,\b-q+J2^k,q,Xcl^qCk- (5.8) 

A; 

The set of renormalization equations has to be solved 
self-consistently: One chooses some values for the expec- 
tation values. With these values, the numerical evalua- 
tion starts from the cutoff A of the original model H and 



proceeds step by step to A = 0. For A = 0, the Hamil- 
tonian and the one-particle operators are fully renormal- 
ized. The case A = allows the re-calculation of all 
expectation values, and the renormalization procedure 
starts again with the improved expectation values by re- 
ducing again the cutoff from A to A = 0. After a sufh- 
cient number of such cycles, the expectation values are 
converged and the renormalization equations are solved 
self-consistently. Thus, we finally obtain an effectively 
free model. 



E, 



(5.9) 



where we have introduced the renormalized dispersion 
relations ik = limA^o £fc,A and ojq — lim^^o'^g.A, and 
the energy shift E = lim^^o Ex ■ 

For the numerical evaluation of the renormalization 
equations we choose a lattice size of = 1000 sites. The 
temperature is fixed to T = 0. 



B. Adiabatic case 

At first, let us discuss our results for the so-called adia- 
batic case Wo «C t. They are shown in panel (a) of Figs. O 
miSl and in panels (a) and (b) of Fig. [6l First, according 
to Fig. [3^ the phononic quasi-particle energies uiq (half- 
filling) are found to gain dispersion due to the coupling 
between electronic and phononic degrees of freedom in 
particular around q = it. Furthermore, if the coupling 
exceeds a critical value gc non-physical negative energies 
at q = TT occur. This feature signals the break-down 
of the present description for the metallic phase at the 
quantum-phase transition to the insulating Peierls state. 

Whereas at half-filling the phonon softening occurs 
at the Brillouin-zone boundary, soft phonon modes are 
found at 2kF = 27r/3 and at 2kF = tt/2 for filling 1/3 and 
1/4, respectively. This can be seen in Fig. [6l Since the 
phonon softening can be considered as a precursor effect 
of the metal-insulator transition, the type of the broken 
symmetry in the insulating phase strongly depends on the 
filling of the electronic band. Note that the critical EP 
coupling gc of the phase transition may be determined 
from the vanishing of the phonon mode (see Ref. [2^. At 
half- filling and for ujo = OAt, a value of gc = 0.31t is 
found, which is somewhat larger than the DMRG result 
of gc = 0.28t of Refs. \M a.nd\M In subsection IVLXI the 
determination of the critical coupling gc within our PRM 
approach will be discussed in more detail. 

Fig. 2^ shows the phonon distribution — (6j6g) for 
the same parameter values as in Fig. |3^. There are two 
pronounced maxima found at wave numbers q — t: and 
q « 0. The peak at q = tt is directly connected to the 
softening of cjq at the zone boundary and can therefore be 
considered as a precursor of the transition to a dimerized 
state. For the critical EP coupling g — gc a divergency 
of should appear at g = tt. The second peak around 
q w follows from renormalization contributions which 




FIG. 3 (Color online) Bosonic quasi-particle energies Uq/uoo 
at half-filling as function of q for different values of the EP cou- 
pling g in the adiabatic case tuo/t — 0.05 (panel (a)), the in- 
termediate case ujo/t = 2.8 (panel (b)), and the anti-adiabatic 
case ijo/t = 6.0 (panel (c)). 



FIG. 4 (Color online) Phonon distribution — {h\bq) as 
function of q for the same parameters as in Fig. [S] 



C. Intermediate case 



become strong for small q for the adiabatic case Wf) <^t. 
This will be explained in more detail in the discussion 
part below. 

Finally, in Fig. [5^ the renormalized fermionic one- 
particle energy Sk is shown in relation to the original 
dispersion Sk = — 2icosfca for the same parameter val- 
ues as in Fig. [3^. Though the absolute changes are quite 
small, the difference between Sk and Ek is strongest in the 
vicinity of A: = and k = tt. In particular, we find ik < Sk 
for fc = and Sk > Ek for fc = tt, so that the renormal- 
ized bandwidth becomes larger than 4t, i.e. larger than 
the original bandwidth. 



Next, let us discuss the results for phonon frequen- 
cies Wo of the order of the hopping matrix element t (in- 
termediate case). The results are found in the panels 
(b) of Figs. [3l m [5l In contrast to the adiabatic case, 
the renormalized phonon energy ujq (Fig. [3]d) now shows 
a noticeable 'kink' at an intermediate wave vector (for 
ujo/t — 2.8). This particular q value, which will be called 
qk in the following strongly depends on the initial phonon 
energy cjo. The appearance of such a 'kink' at gfe < tt 
is a specific feature of the intermediate case. The wave 
number qk is characterized by a strong renormalization 
of the phonon energy in a small g-range around qk , where 
Cjq/ijj{) > 1 for g < gfe and ujq/uio < 1 ioi q > qk holds. 
The origin of these features will be discussed in more 
detail below. 

Similar to w^, also the phonon distribution in 



21 



0.04 




I 

IW 



— I — 
(a) 



■ - ■ g / 1 = 0.06 

g/t = 0.15 

g/t = 0.24 




1.5 



FIG. 5 (Color online) Fermionic quasi-particle energies (ik — 
ek)/t as function of k for the same parameters as in Fig. [S] 
Here Sk is the original electronic dispersion. 



Fig. 2)3 shows a pronounced structure of considerable 
weight around qk- Finally, in Fig. [5Jd the difference of 
the fermionic one-particle energies (e^ — Ek) is shown. 
Again a remarkable structure is found, though the abso- 
lute changes are small for the present (/-values. 



D. Anti-adiabatic case 

Finally, let us discuss the results for the anti-adiabatic 
case Wo » t. In panels (c) of Figs. [31 IH [5] a value of 
ujQ/t = 6.0 was used. As most important feature a stiff- 
ening of the renormalized phonon frequency ioq (Fig. [3J;) 
is found instead of a softening as in the adiabatic case. In 
particular, for large values of the EP coupling no soften- 
ing of the phonon modes is found at g = tt. Moreover, no 
large renormalization contributions occur in any limited 
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FIG. 6 (Color online) (a) Phononic quasi-particle energy 
ubq in unit of ujo of the one-dimensional HM with 500 lattice 
sites for filling 1/3 and different values of the EP coupling g. 
ujo/t = 0.05. (b) Same quantity Cuq/uo for filling 1/4. 



g-space regime which would lead to peak-like structures. 
Instead an overall smooth behavior is found in the entire 
Brillouin zone. 

Also the phonon distribution (Fig. SI;) shows a 
smooth behavior with a maximum at g = tt. The lack of 
strong peak-like structures in q space indicates that there 
is no phonon mode that gives a dominant contribution to 
the renormalization processes. 

If one compares the renormalized electronic bandwidth 
for the anti-adiabatic case (Fig. (5};) with that of the adi- 
abatic case (Fig. [S^), one observes a relatively strong 
reduction of the bandwidth. This indicates the tendency 
to localization in the anti-adiabatic case. It also indicates 
that the metal-insulator transition in the anti-adiabatic 
limit can be understood as the formation of small im- 
mobile polarons with electrons surrounded by clouds of 
phonon excitations. In the present PRM approach, a 
renormalized one-particle excitation like Sk corresponds 
to a quasiparticle of the coupled many-particle system. 
Therefore, a completely flat k dependence of would be 
expected to be found in the insulating regime. 



E. Discussion 

It may be worthwhile to demonstrate that the PRM 
approach has the advantage that all features of the results 
for Li)q and or can easily be understood on the basis 
of the former renormalization equations. For simplicity. 
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we shall restrict ourselves to the case of half-filling and 
to the renormalization of the phonon energies (bq. 

The basic equation is the renormalization equation 
(|5.6p . Due to the 0-functions 8fc^5(A, AA) in all equa- 
tions a renormalization approximately occurs when the 
energy difference |wg,A + £fc,A — Sk+q,\ \ lies within a small 
energy shell between A and A — AA. As one can see from 
(15. 6p the most dominant renormalization processes take 
place for small values of the cutoff A. Therefore, the 
largest renormalization contributions come from fc and q 
values that fulfill the condition 

efc+g,A - Efc^A ~ Wg,A- (5.10) 

From (|5.6p directly follows a second condition for the 
renormalization contributions to LOq.\. Due to the expec- 
tation values (n\ — n'f,^^) in (|5.6p the renormalization of 
ujq^x is caused from the coupling to particle-hole excita- 
tions. Therefore, the energies Sk^x and Sk+q.x have to be 
either below or above the Fermi level, i.e. |fc| < kp and 
|A; -I- gl > kp or \k\ > kp and |fc -I- q| < kp. 

Let us first discuss the adiabatic case wq <C t. The 
most dominant contributions to the renormalization are 
expected when both conditions are simultaneously ful- 
filled. This is the case for q w ±7r or partially also for 
q « 0. Note that for q ~ tt practically all fc-values 
can contribute to the renormalization of (j5.6p . which 
is not the case for q-values different from tt. For in- 
stance, for q only few k points from the sum in 
(15. 6p can contribute which are located in a small region 
around the Fermi momentum kp. On the other hand, for 
g « 0, the energy denominator is almost zero so that still 
some noticeable renormalization structures are found in 
Fig. [3^. Moreover, for the adiabatic case, where u!q^\ is 
small, the energy denominator of (j5.6p can be replaced 
by (£fc,A — Sk+q,\)- Therefore, almost all particle-hole 
contributions to w^.a are negative because {n^. — n'f.^^) 
and (efe,A — £k+q,x) have always different signs. One con- 
cludes that in the adiabatic case ujq^x will be renormal- 
ized to smaller values where the renormalization at q — tt 
should be dominant. 

The behavior of Qq for the case of intermediate phonon 
frequencies {ujo/t = 2.8 in Fig. [3]3 and Fig. [Dd) can again 
be understood on the basis of the renormalization equa- 
tions (|5.6p and condition (|5.10p . As was already dis- 
cussed, particle-hole excitations lead to the renormaliza- 
tion of w<j,A- Therefore, from the sum over k in Eq. (|5.6p 
only k terms contribute where either |fc| < kp and 
|fc + g| > kp or |fc| > kp and |fc -f qj < kp. For the latter 
case always (e^^A — £k+q,x) > is vahd so that (|5.10p 
can not be fulfilled. Therefore, we can restrict ourselves 
to contributions |fc| < kp und |fc + g| > kp, for which 
always {ek,\ - ek+q,x) < and {nl - nl^^) > holds. 
The largest renormalization should result from a small q 
region around some q vector qk for which Sk+q^ ~£k = ^j->a 
is approximately fulfilled. Since is of the order of t, 
qk is located somewhere in the middle of the Brillouin 
zone and depends strongly on loq. From Eq. (|5.6p also 
follows that renormalization contributions to ujq change 



their sign at qk due to the sign change in the energy de- 
nominator. 

Finally, from equation (j5.6p one may point out also the 
stiffening of the phonon modes in the anti-adiabatic case 
ujo/t = 6.0. In this case the phonon energy loq is much 
larger than the electronic bandwidth. Therefore, for all 
A a positive energy denominator {ujq,\ + Sk.x — £fc+<;,A) is 
obtained. Nevertheless, for half-filling in the k sum on 
the right hand side of (j5.6p there are as many negative 
as positive terms due to the factor (n^ — ^iJc+g)- Since 
from {nl - nl^^) < always follows (e^^A - £fc+g,A) > 0, 
the negative terms have larger energy denominators and 
are always smaller than the positive terms. The result- 
ing renormalization of Wq^\ is therefore positive for all q 
values and largest for q = n due to the smallest energy 
denominator. 



VI. QUANTUM PHASE TRANSITION IN THE 
ONE-DIMENSIONAL HOLSTEIN MODEL 

In this section we want to demonstrate the ability of 
the PRM approach to describe also quantum phase tran- 
sitions. In particular, we shall investigate the transition 
from the metallic to the insulating charge ordered phase 
when the electron-phonon coupling g exceeds a critical 
value. 



A. Uniform description of metallic and insulating phases 
at half-filling 

In the following we present a uniform description that 
covers the metallic as well as the insulating phase of the 
HM in the adiabatic case. We mainly follow the approach 
of Ref. m where we have discussed methodological as- 
pects in more detail. As already mentioned above, the 
simple approach of subsection IV. Al breaks down for EP 
couplings g larger than some critical value gc where a 
long-range charge density wave occurs and the ions are 
shifted away from their symmetric positions. An ade- 
quate theoretical description needs to take into account 
a broken symmetry field. For this purpose, the under- 
lying idea of subsection IIII.CI to take such a term into 
account in the renormalization ansatz will be transferred 
to the present case. As one can see from Fig. [51 the or- 
der parameter of the insulating phase strongly depends 
on the filling of the electronic band. Therefore, in the 
following we restrict ourselves to the case of half-filling. 
Here, the unit cell is doubled and a dimerization occurs 
in the insulating phase. 

Following Ref. [s^, the Hamiltonian in the reduced Bril- 
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louin zone including symmetry breaking fields reads 
Hx = Tio^x+Hi^x, (6.1) 

k>0,a (J>0,7 

+ Ex + Y, ^fc,A (cS,fcCi,fe + h.c.) 

k 

^ A:,?>0 

where AJ: ^ and are the appropriate order parameters 
for the electronic and the phononic symmetry breaking 
fields. Note that the reduced Brillouin zone leads to ad- 
ditional band indices a, /?, 7 = 0, 1 of both electronic and 
phononic one-particle operators. Furthermore, we de- 
fined 6 A = A ~ (A) and Q = w/a. The ansatz (|6.ip is 
restricted to the one-dimensional case at half-filling. To 
extend the approach to higher dimensions one would need 
to take into account all Q wave vectors of the Brillouin 
zone boundary. 

Before we can proceed we need to diagonalize 7io,A- 
For this purpose a rotation in the fermionic subspace 
and a translation to new ionic equilibrium positions are 
performed in order to diagonalize 7io,A 

T-Co.x ^ ^^<^a,k,xCi^k,x'^a,k,x (6.2) 

k>0 a 
q>0 7 

with new fermionic and bosonic creation an annihila- 
tion operators, cj^\ x ^^"^ ^7^9 A' rewrite Tii^x 
in terms of the new operators, c'^\ ^ and B^]^ ^. 

Finally, we have to transform Hx to ?i (a-aa) according 
to (|2.17p to derive the renormalization equations for the 
parameters of Ti.x. Here the ansatz 

k,q a.f^jj 

x{SBlJ{Cl^Cf,,k+,,x)~h.c.] 

is used. The coefficients ^^'^^x aa ^^'^^ to be fixed in such 
a way so that only excitations with energies smaller than 
(A — AA) contribute to Ti-ij^x-AX)- The renormalization 
equations for the parameters ea,k.x, A^ ^,ujy,q^x, A^, and 
dk'q 'x finally obtained by comparison with (|6.ip af- 
ter the creation and annihilation operators C^^\ S^^^ ^ 
have been transformed back to the original operators 
cl}\,bj}q- The actual calculations are done in close anal- 
ogy to subsection IV. Al Note that again a factorization 
approximation was used and only operators of the same 



structure as in (|6.ip are kept. Therefore, the final renor- 
malization equations still depend on unknown expecta- 
tion values, which are evaluated with the full Hamilto- 
nian Ti.. Note that in order to evaluate the expectation 
values (A) = {Ax)hx additional renormalization equa- 
tions have also to be found for the fermionic and bosonic 
one-particle operators, c]^ and bt^ ^. By using the same 
approximations as for the Hamiltonian a resulting set of 
renormalization equations is derived. It is solved numer- 
ically where the equations for the expectation values are 
taken into account in a self-consistency loop. 

By eliminating all excitations in steps AA we finally 
arrive at cutoff A = which again provides an effectively 
free model H — limA^o Hx = liniA^o Wo, a. It reads 

^ = H ,^Ca,fc + ^ A^ (cS,fcCi,fc + h.c.ye.S) 

k>0,a k>0 
9>0,7 

where it was defined Sa.k = liniA^o £Q,fe,A7 A| = 
liniA-^o A^^^, Lo^^q = limA^o ^^7,9, a, and A'' = limA^o A^;. 
Note that all excitations from TLi^x were used up to 
renormalize the parameters of TLq. The expectation val- 
ues are also calculated in the limit A — > 0. Because 
7i is a free model they can easily be determined from 
{A)-H = {Ax)h>. = ((IhnA^o Ax)).H ■ 

B. Results 

In the following, we first demonstrate that the PRM 
can be used to investigate the Peierls transition of the 
one-dimensional spinless Holstein model (|5.ip at half- 
filling. The phonon energy is fixed to = O.lt. In 
particular, our analytical approach provides a simultane- 
ous theoretical description for both the metallic and the 
insulating phase. Finally, we compare our results with 
recent DMRG calculations [H; HI. 

First, let us consider the critical electron- phonon cou- 
pling (7c- For that purpose, in Fig. [7] a characteristic elec- 
tronic excitation gap A for infinite system size is plotted 
as function of the EP coupling g, where A was deter- 
mined from the opening of a gap in the quasi-particle 
energy Sk (see text below). A closer inspection of the 
data shows that an insulating phase with a finite excita- 
tion gap is obtained for g values larger than the critical 
EP coupling gc « 0.24t. A comparison with the criti- 
cal value gc « 0.28t obtained from DMRG calculations 
I43I : I45I shows that the critical values from the PRM ap- 
proach might be somewhat too small. However, this dif- 
ference can be attributed to the exploited factorization 
approximation in the PRM which suppresses fluctuations 
so that the ordered insulating phase is stabilized. Note 
that in order to determine gc a careful finite-size scaling 
was performed as shown for some g values in the inset 
of Fig. [71 A linear regression was applied to extrapolate 
our results to infinite system size. Note that the finite 
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FIG. 7 Electronic excitation gap of the one-dimensional HM 
at half-filling where the data are extrapolated to an infinite 
chain. The inset shows the finite-size scaling for g values of 
the EP coupling of 0.26* (circles), 0.266t (triangles), and 0.27t 
(diamonds). 



size scaling may be affected by two different effects: Sup- 
pression of long-range fluctuations by the finite cluster 
size and by the used factorization approximation so that 
a rather unusual dependence on the system size is found. 

In contrast to other methods, the PRM directly pro- 
vides the quasi-particle energies: After the renormahza- 
tion equations were solved self-consistently the electronic 
and phononic quasi-particle energies of the system, Sk 
and ujq, respectively, are given by the limit A — > of the 
parameters ej^ j, and q \ of the diagonal Hamiltonian 
7^0, (A^o) of (|6.2p . In Fig.[S]the renormalized one-particle 

energies ik = £a=o,k,x=o and ujq = w^^o,g,A=o as quasi- 
particle of the full system are shown for different values 
of the EP coupling g. The upper panel shows that the 
electronic one-particle energies depend only slightly on g 
as long as g is smaller than the critical value gc ~ 0.24i. 
If the EP coupling g is further increased a gap A opens 
at the Fermi energy so that the system becomes an insu- 
lator. Remember that the gap A has been used as order 
parameter to determine the critical EP coupling g^ of the 
metal-insulator transition (see Fig. [7]). The lower panel of 
Fig. [8] shows the results for the phononic one-particle en- 
ergy ibq. One can see that Qq gains dispersion due to the 
coupling g between the electronic and phononic degrees of 
freedom. In particular, the phonon mode at momentum 
2fci?, i.e. at the Brillouin-zone boundary becomes soft if 
the EP coupling is increased up to gc ~ 0.24t. However, 
in contrast to the metallic solution of subsection IV. Al Cj„ 
at 2kF always remains positive though it is very small. 
Note that for g values larger than g^ the energy Qq in- 
creases again. This phonon softening at the phase tran- 
sition has to be interpreted as a lattice instability which 
leads to the formation of the insulating Peierls state for 




q/7l 

FIG. 8 (Color online) Fermionic quasi-particle energy ik = 
£ci=o,fc A=o (upper panel) and bosonic quasi-particle energy 
ujq — uj^^o q x^o (lower panel) of a chain with 500 lattice sites 
for different EP couplings g. 



g > gc- The phase transition is associated with a shift 
of the ionic equilibrium positions. A lattice stiffening oc- 
curs if g is further increased to values much larger than 
the critical value gc ~ 0.24i. 

Note also that the critical coupling gc ~ 0.24t ob- 
tained from the opening of the gap in ik is significantly 
smaller than the gc value of ~ 0.31t which was found from 
the vanishing of the phonon mode at the Brillouin zone 
boundary in the metallic solution of subsection IV. Al In- 
stead, one would expect that both the gap in and the 
vanishing of ujq should occur at the same gc value. This 
inconsistency can again be understood from the factor- 
ization approximation in the PRM: As discussed above, 
the inclusion of additional fluctuations leads to a less sta- 
ble insulating phase so that a gc value larger than 0.24t 
would follow. On the other hand, the dispersion of ujq 
due to renormalization processes would be enhanced by 
taking additional fluctuations into account. Thus, a gc 
value smaller than « 0.31i would follow. In this way, 
both ways to determine gc would be consistent with each 
other and could lead to a common result for gc in be- 
tween 0.24t and 0.31t. This would be in agreement with 
the DMRG value of gc ~ 0.28t (jH; SI). 
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VII. CHARGE ORDERING AND SUPERCONDUCTIVITY 
IN THE TWO-DIMENSIONAL HOLSTEIN MODEL 

As a second example for a quantum phase transi- 
tion, we now study the competition of charge-density 
waves (CDW) and superconductivity (SC) for the 
two-dimensional half-filled Holstein model by use of 
the projector-based renormalization method. In one 
dimension the coupling of electrons to phonons gives 
rise to a metal-insulator transition. In two dimensions 
the electron-phonon interaction may also be responsible 
for the formation of Cooper pairs. In the following, the 
competing influence of superconductivity and charge or- 
der will be discussed for two dimensions. The PRM not 
only allows to study SC and CDW correlation functions 
but gives direct access to the order parameters. The 
discussion closely follows the approach of Ref. [H 

The relationship between a possible superconducting 
and an insulating Peierls-CDW phase in the 2d-Holstein 
model has been subject to a number of studies in the 
literature (for details we refer to Ref. H^). In general, it 
is believed that the onset of strong SC correlations sup- 
presses the development of CDW correlations and vice 
versa. Thus close to the phase transition, both types of 
correlations must be taken into account. 



A. Unified description of SC and CDW phases at half-filling 

To find a uniform description of both the supercon- 
ducting (SC) and the insulating CDW phase, two fields, 
which break the translation and the gauge symmetry 
should be added to the Hamiltonian. Thus, the model 
on a square Lattice is given by 

n = Hq+Hi (7.1) 

k,cr q 

+ E(^k4,T^-k,i + Arc-kaCk,T) 

k 

+ ^ E (^k 4,.Ck-Q,. + h.c) + VNA\b^^ + 6q) 

k,CT 

^1 = E {^H,^Ck+q.<T+ Vk+q.<xCk.<T}. (7.3) 

k.q.cT 

where k is the wave vector on the reciprocal lattice and 
Q is the characteristic wave vector of the CDW phase 
Q = (7r/a,7r/a). Assuming an electron hopping between 
nearest-neighbor sites, the electronic dispersion is given 
by £k = — 2t(cos kxtt + cos kyu) — /i, where /i is the chem- 
ical potential. Moreover, ujq is the dispcrsionless phonon 
energy, and g denotes the coupling strength between the 
electrons and phonons. At the beginning of the renormal- 
ization the two symmetry breaking fields and Aj^, 



as well as A*", are assumed to be infinitesimally small 
(A^^O, AP^CA'-^ O). ^ 

The unperturbed Hamiltonian Tio can be diagonalized, 
since its electronic part is quadratic in the fcrmionic op- 
erators. Note that due to the doubling of the unit cell 
in the insulating phase, in TYo the creation operator cj^. ^ 
is coupled to Ck-q.o-- In addition the coupling of cj^ ^ 

to c^jj. I is caused by superconductivity. Therefore, the 
eigenmodes of Tio can be represented as a linear combi- 
nation of the following four operators 

( c-k-qa 4,T c-k,i 4+Q,| ) (7.4) 

In the renormalization procedure, all transitions with 
energies larger than A will be integrated out. As can be 
seen, the renormalized Hamiltonian can again be divided 
into Hx = Ho, A + 'Hi, A- If one denotes hy aj^j^-^ {ce — 
1 • • • 4) the A dependent eigenmodes of 7io,A the electronic 
part of the renormalized Hamiltonian Tio.Acan be written 
as 

Ho'a E {-^l.k.A (aI,k.Aai,k,A + a2,k,A«2,k,A) 

ker.BZ 

+£^2,k,A (a3,k,A"3,k,A + a4_k,A"4,k,A^ | (7.5) 

where the eigenenergies are given by 

i^i/2,k,A - '•--^ ± M^k,A (7.6) 

^^k,A = ^( ''^'^'f-^-^ )' + |ALAP + |ALAP 

for ek,A + £k-Q,A > 0, whereas for £k,A + £k-Q,A < the 
±-signs have to be reversed. Note that in ()7.6p the sum 
of the two order parameters squared enter the energies 

Ea,k,X of CM- 

In order to derive the renormalization equations, the 
unitary transformation (j2.17p has to be evaluated explic- 
itly. Thereby, also the interaction Hi^\ has to be ex- 
pressed in terms of the eigenmodes aA,k,A of Ho, a- More- 
over, an ansatz for X^ /^x has to be made in analogy 
to what was done in the previous sections. The explicit 
calculation is found in Ref. [s^ . 

B. Results and Discussion 

For the numerical evaluation of the renormalization 
equations, we consider a square lattice with N — 144 
sites. The temperature is set equal to T = 0, and a small 
value of LOo = O.li is chosen. For simplicity, we also re- 
strict ourselves to s-wave-like superconducting solutions. 

The results are shown in Fig. [51 where the k-dependent 
symmetry breaking fields A^ (black) and A^ (red) for 
k — (7r/2,7r/2) are plotted as function of the electron- 
phonon coupling g. The coupling g is restricted to small 
values g/2t < 0.04. As can be seen from Fig. [51 for small 
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0.04 



g/2t 

FIG. 9 (Color online) Renormalized values of the Peierls gap 
Aj^ (black line) and of the superconducting gap A|. (red line) 
at wave vector k = (7r/2,7r/2). A square lattice with 144 
lattice sites at half-filling was taken, tJo/f = 0.1 and T — 0. 



values of g/2t < 0.010 the system is in a pure supercon- 
ducting state, i.e. no charge order is present. For small g, 
the superconducting gap increases roughly proportional 
to g^. In the intermediate g range, 0.010 < g/2t < 0.023, 
a coexistence of both phases is found. The system is in a 
combined superconducting-charge ordered phase. Here, 
the g dependence of is no longer quadratic as in the 
small g regime. Instead, reaches a maximum value 
and drops down to zero with increasing g. Finally, for 
g/2t > 0.023 the superconducting phase is completely 
suppressed and the system is in a pure charge ordered 
state. 



VIII. SUMMARY 

The aim of this contribution was to discuss the ba- 
sic ideas of a new theoretical approach for many-particle 
systems which is called projector-based renormalization 
method (PRM) and its application to a number of non- 
trivial physical problems. Instead of eliminating high- 
energy states as in usual renormalization group meth- 
ods in the PRM high-energy transitions are successively 
eliminated. Thereby, a unitary transformation is used 
where all states of the unitary space of the interacting 
system are kept. In that respect, the PRM is closely re- 
lated to the similarity transformation introduced by Wil- 
son and Glazek and to Wegner's flow equation method 
though both approaches start from a continuous formula- 
tion of the unitary transformation. The PRM starts from 
a Hamiltonian which can be decomposed into a solvable 
unperturbed part and a perturbation, H — Ho + Hi, 
where the latter part induces transitions between the 
eigenstates of Ho . 

Suppose a renormalized Hamiltonian H\ has been con- 
structed which only contains transitions with transition 



energies smaller than some given cutoff energy A. The 
Hamiltonian H\ can be further renormalized by elimi- 
nating all transitions from, roughly speaking, the energy 
shell between the cutoff A and a reduced cutoff (A — AA), 
and so on. This is done by a unitary transformation 
H(x-A\) — e^^'^^ H\ e'-^^'^^ which guarantees that the 
eigenspectrum is not changed. The generator of the uni- 
tary transformation X^.aa is specified by the condition 
Qa-aa'^(a-aa) = where Qa-aa is the projector on all 
transitions with energy differences larger than (A — AA). 
The latter condition implies that all transitions from the 
'shell' between A and A — AA are eliminated and lead 
to a renormalization of 7i(A-AA)- Note that only the 
equivalent part Qa-aa-^a,aa of Xa,aa is fixed whereas 
the orthogonal part Pa-aa-'^'a.aa can be chosen arbi- 
trarily. Note that this additional freedom can be used in 
a different way. Whereas in the original version of the 
PRM the remaining part Pa-aa-'^a.aa of X\^ax was set 
equal to zero for simplicity this part was used in Weg- 
ner's flow equation method as the only relevant part when 
the transformation was performed continuously. In this 
case, the interaction parameters were chosen to decay ex- 
ponentially. By proceeding the renormalization up to the 
final cutoff A = all transitions induced by Hi.x are elim- 
inated. The final renormalized Hamiltonian H — Ho.\=o 
is diagonal and allows to evaluate in principle any cor- 
relation function of physical interest. In particular the 
one-particle excitations of H can be considered as quasi- 
particles of the coupled many-particle system since the 
eigenspectrum of the original interacting Hamiltonian H 
and of H are in principle the same since both are con- 
nected by a unitary transformation. 

Note that the present approach has the advantage 
of formulating the renormalization quite universally. 
By specifying the unitary transformation of the many- 
particle system both the PRM and Wegner's flow equa- 
tion method can be derived from the same basic ideas. 
However, the stepwise transformation of the PRM has 
its own merits. Firstly, as was shown in Sec. Illl.Ci 
Sec. lVIl and Sec. I VIII the physical behavior on both sides 
of a quantum critical point can be described within the 
same PRM scheme. This seems not the case for the flow 
equation approach. In particular, by allowing symme- 
try breaking terms in the 'unperturbed' part Hq^x, the 
transformation of eigenmodes of the Liouville operator 
Lo,A can be followed in each renormalization step. This 
makes the description of quantum critical points possi- 
ble. Secondly, in Sec. III.BI a perturbation theory for Hx 
was given. This allows to evaluate physical properties in 
perturbation theory. In contrast to a recent perturbation 
approach on the basis of the flow equation method, in the 
PRM no equidistant spectrum of Ho is required. 
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TABLE I Action of the 7^(j) as used in the calculations. For 
convenience, the dimer indices of the states are suppressed. 



APPENDIX A: Example: dimerlzed and frustrated spin 
chain 

In this appendix we are going to investigate ground- 
state properties of a dimerized and frustrated spin chain. 
We apply the projector-based perturbation theory and 
use expression (j2.15|) for Hx and chose A ^ right from 
the beginning. In this case, the interaction Tii is com- 
pletely integrated out in one step. The starting Hamil- 
tonian reads 

n = Ho+Hi, (Al) 

^0 — J S2iS2i+l, 

i 

'Hi = J y^[aS2iS2i-i + /3 (S2tS2t-2 + 82^-1821+1)], 





\t^.t^) 




-J{a-2P)\s,f°-^) 
J(a + 2/3) 

./(a + 2/5){|t»,t»)-|t±,iT)} 
J(a + 2P) {\t+,t-) + \t-.t+)} 




\s,t^)..r,s) 
\s,f),\t°,s) 
\s.t-),\t-,s) 




Ja {\t+,fi) - \fi..t^)} 
Ja {\t+,t-) - \t-.t+)} 
Ja {\t°,t-) - \t-.t")} 




\s,s) 




J {a - 2/3) {\t+J.-) + \f-,t+) 



The model itself is of some physical interest because it 
can be used to describe some spin-Peierls compounds like 
CuGeOa or TTFCuBDT 

In the following we are interested in the limit of strong 
dimerization of the model so that we start from isolated 
dimers as described by Hq. Every dimer can be in the 
singlet state or in one of the three degenerated triplet 
states. Here, the dimer states are energetically separated 
by the singlet-triplet splitting, A = St — Sg — J- Thus, 
triplets can be considered as the basic excitations of the 
system. 

Following the ideas of Refs. [3 and [T9I, the contribu- 
tions to the perturbation Tii can be classified according 
to the number of created or annihilated local triplets. 

Hi = (A2) 

j 

The introduced excitation operators %nij) only act on 
the local dimers with indices j and j — I, and create 
TO local triplets. The %n{j) are eigenoperators of the 
Liouville operator Lq and the corresponding eigenvalues 
are — mA. The actual contributions to To^j), Ti{j), 
and T2ij) are summarized in Table IH and T-i{j) and 
T_2(j) are given by the relation T_m{j) = [Tm{j)]'^- 

In the limit of strong dimerization, Hilbert space sec- 
tors with different numbers of triplets in the system are 
energetically separated because the unperturbed part TLq 
of the Hamiltonian (|Aip does not change the number 
of triplets in the system, and the interaction TLi only 
leads to modest corrections. Consequently, the evalua- 
tion of the effective Hamiltonian (|2.15p can be simplified 
if one concentrates on a Hilbert space sector with a given 
fixed number of triplets. In the following, actual cal- 
culations are presented for the two energetically lowest 
sectors where the system contains no or only one triplet. 



The subspace without triplets consists of a single state, 
i.e. the singlet product state, |<i>Gs) = 1*2) ■ • ■ \sn)- 



Because the effective Hamiltonian Ti. 



(A-.0) 



is obtained 



from the original Hamiltonian Ti. by means of a unitary 
transformation, the ground-state energy can be calcu- 
lated from 



E'gs = lini {Ti) = hm 

13 — 'oo /3 — ^(x 



Tr?^ 



(A-»0) 



Tr e~''^(^^o) 



(A^O) I'I'GS 



Here U 



(A-,0) 



is given by (|2.15p where P(;y^o)^i ' 
J2j'^o{j) and Q(_)^^g)Hl is the remaining part of (|A2 
Using the notation of Ref. [l^ one easily finds 



GS 



= NJ 



3 3 , 
(a 

4 32^ 



0{nl). (A3) 



This result agrees with findings of Refs. [ij andll9|. Note 
that higher order terms can easily be calculated by imple- 
menting a computer based evaluation algorithm as dis- 
cussed in Ref. 19^ where a cumulant method (fisl ) was 
applied to the same model. 

The case of a single triplet in the system is more com- 
plex because a triplet can easily move along the chain. 
Consequently, it is advantageous to introduce momentum 
dependent states, 



\sn) 



and the eigenvalues of this Hilbert space sector can be 
calculated by = limA^o {^k\ 'Hx \^k)- We again em- 
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ploy the useful notation of Ref. and obtain 



El 



1 



1 



(A4) 



J I ^ (a - 2/3) + I cos(fca) 

- — J (a - 2pf cos(2fca) + 0{HI). 
16 

Note that the energy gap of the system can easily be de- 
termined from Eq. (|A4p by considering the case fc = 0. 
The k dependence of Ej^ describes the triplet dispersion 
relation. Furthermore, the calculation can easily be ex- 
tended to higher orders. 

The same model was also studied based on Weg- 
ner's flow equation method (^) where both ground-state 
energy and triplet dispersion relation were calculated in 
high orders. However, for this purpose a set of coupled 
differential equations had to be integrated so that this ap- 
proach is restricted to systems with an equidistant eigen- 
value spectrum of the unperturbed part Tip of the Hamil- 
tonian. 
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